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SUMMARY 


A computer program has been developed which is capable of calculating the 
flow field in the supersonic portion of a mixed-compression aircraft inlet 
operating at angle of attack. The calculation procedure is based on the 
method of characteristics for steady three-dimensional flow. The bow shock 
wave and the internal shock' wave system are computed using a discrete shock 
wave fitting procedure. 

The computer program has the capability to compute the internal flow 
field with or without the discrete fitting of the internal shock wave system. 

The option in which the internal shock wave system is not fitted can be em- 
ployed in situations in which the strength of the internal shock wave system 
is weak,. and thereby an acceptable solution can be obtained by smearing the • 
internal discontinuities. 

The influence of molecular transport can be included in the computation 
of the external flow about the forebody, and in the computation of the inter- 
nal flow in which the shock waves are not discretely fitted. This is accom- 
plished by treating the viscous and thermal diffusion terms in the governing 
partial differential equations as forcing functions, or correction terms, in 
the method of characteristics scheme. 

The thermodynaniic model employed in the computer program is contained in 
a separate group of subroutines. The assumed thermodynamic model is that of a 
thermally and calorically perfect gas. Other thermodynamic models may be 
employed by suitably modifying the existing subroutines or by replacing them. 

The representations for the molecular transport properties, viscosity and 
thermal conductivity, are contained in a separate subroutine. Dynamic viscosity 
is represented by Sutherland's law. Thermal conductivity is represented as a 



quadratic function of temperature. Alternative formulations for the transport 
properties, may be employed by modifying the existing subroutine or by replacing 
it. 


The contours of the centerbody and the cowl are represented by a separate 
subroutine. The existing subroutine has the capability to describe a variety 
of axisymmetric contours. Other geometries, such as those having elliptic or 
Super-elliptic cross-sections, may be described by suitably modifying the 
existing subroutine or by -replacing it. 

A major assumption of the present analysis is that the cowl Tip is con- 
tained in a plane of constant x. Moreover, it is assumed that both the center- 
body contour and the cowl contour are smooth and have continuous first partial 
derivatives. 

The computer program cannot: 

1. compute subsonic flow, 

2. compute the external flow field about the forebody if the 

bow shock wave does not exist entirely around the forebody, or 

3. compute the internal flow field if the bow shock wave is ingested 
into the arinulus 

The computer program is written in Fortran IV for the CDC 6500 computer. 
The program can be easily modified to be compatible with other computers. 
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SECTION I 


INTRODUCTION 


A computer program has been developed for calculating the flow field in 
the supersonic portion of a mixed-compression aircraft inlet operating at angle 
of attack. The general features of the inlet geometry and the flow field are 
illustrated in Figure 1. The theoretical analysis on which the computer pro-' 
gram is based is presented by Vadyak, Hoffman, and Bishop (1). 

This report presents a discussion of the computer program organization, de 
scriptions of the subroutines, a discussion of the input parameters, a brief 
interpreta tion of the output information, and five sample cases to illustrate 
the application of the computer program. 

A detailed description of the theoretical analysis used in this program is 
given in NASA CR- 135425, "The Flow Field Calculation in Supersonic Mixed- 
Compression Aircraft' Inlets at Angle of Attack Using the Three-Dimensional 
Method of Characteristics with Discrete Shock Wave Fitting" by Joseph Vadyak 
and Joe D. Hoffman. 
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EXTERNAL FLOW REGIME 


INTERNAL FLOW REGIME 



FIGURE L MIXED-COMPRESSION AIRCRAFT INLET 



SECTION II 


PROGRAM ORGANIZATION 


1 . INTRODUCTION 

In this section, the program overlay structure is presented, and some 
general comments concerning program input and parameter initialization are 
made. The subroutines are identified which contain the thermodynamic model, 
the transport property representations , and the formulations used for repre-’ 
senting the centerbody and the cowl contours. The flow field integration 
options are briefly discussed. Some comments concerning program output are 
presented. 


2. INTEGRATION OPTIONS 

The computational flow regime is divided into two subregimes: the exter- 
nal flow field and the internal flow field (see Figure 1). If desired, only 
the external flow field may be computed. Alternatively, if the solution is 
known at the axial station of the cowl lip, the internal flow field may be 
determined without employing the forebody flow integration option. 

The program has the capability to compute the internal flow field with 
or without the discrete fitting of the internal shock wave system. The option 
in which shock waves are not discretely fitted might be employed if the 
strength of the internal shock wave system is relatively weak, and thereby an 
acceptable solution can be obtained by smearing all internal discontinuities. 

The analysis includes the influence of molecular transport in the com- 
putation by treating the viscous and thermal diffusion terms in the governing 
partial differential equations as forcing functions, or correction terms, in 
themethod of characteristics scheme. At present, the program has the capa- 
bility to include the influence of molecular transport in the computation of 
the external flow field about the forebody, and in the computation of the 
internal flow field in which the shock waves are not discretely fitted. The 
program option in which the internal shock waves are discretely fitted does 
not have the capability to include the influence of molecular transport in 
the computation, but rather assumes the flow to be in viscid and adiabatic. 


3. PROGRAM OVERLAY STRUCTURE 

The overall program consists of .61 routines (6 program routines and 55 
subroutines). The program is too large to be stored continuously in computer 
memory, hence, an overlay scheme is employed. The overlay structure is pre- 
sented in Figure 2. Three levels of -overlaying are used: the resident 
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overlay level, the primary overlay level, and the secondary overlay level. 
0VERLAY (0,0) is the resident overlay and controls the overall execution of the 
program. This overlay contains many of the subroutines which describe the 
thermodynamic model, the molecular transport properties, and the centerbody 
and the cowl contours. 0VERLAY (0,0) calls the two primary level overlays, 
0VERLAY (1,0) and 0VERLAY (2,0). 0VERLAY (1,0) performs data input, parameter 
initialization, and, if desired, internal generation of the initial-value 
plane. 0VERLAY (2,0) is used in conjunction with the secondary overlays, 
0VERLAY (2,1), 0VERLAY (2,2), and 0VERLAY (2,3), to perform the actual flow 
field integration. The primary level overlay 0VERLAY (2,0) contains the 
preponderance of the unit processes and the interpolation schemes. 0VERLAY 
(2,1), which is a secondary level overlay, contains the control logic used in 
performing the integration for the external flow field about the forebody. 

The secondary level overlay 0VERLAY (2,2) contains the control logic used in 
determining the internal flow field in which shock waves are discretely fitted. 
The secondary level overlay 0VERLAY (2,3) contains the control logic used in 
determining the internal flow field in which shock waves are not discretely 
fitted. 

Routine LINKOO, the main control routine in 0VERLAY (0,0), first calls 
0VERLAY (1,0) for data input, parameter initialization, and, if desired, 
internal generation of the initial-value plane. Control is then returned to 
the resident overlay, 0VERLAY (0,0), 0VERLAY (2,0) is then called from LINKOO. 
The main control routine in 0VERLAY (2,0) is routine LINK20. If the forebody 
flow field is to be computed, LINK20 calls 0VERLAY (2,1). Then, if the inter- 
nal flow field is to be computed with discrete fitting of the internal, shock- 
wave system, 0VERLAY (2,2) is called. If the internal flow field is to be 
computed without the discrete fitting of the internal shock wave system, 
0VERLAY (2,3) is called. It should be noted that 0VERIAY (2,2) and 0VERLAY- 
(2,3) both cannot be used in the same computer run. 

The computer program has a restart capability in which an internally 
generated restart file is retrieved from tape. When the restart option is 
used, control is returned to the integration control secondary level overlay 
which was in use at the time the initial execution was terminated [e.g., if 
the forebody flow was being computed, control is returned to 0VERLAY (2,1)]. 


4. DATA INPUT 

The input data are entered through both namelist input and formatted read 
statements. Seven namelists must always be specified: namelists LIST! to 

LIST7. Namelist LIST1 contains the input parameters which specify the flow 
field integration options, the number of planes of symmetry, the output options, 
etc. Namelist LIST2 contains the parameters which specify the free-stream 
conditions, the vehicle orientation, and parameters associated with the speci- 
fication of the initial -value plane. -Namelist LIST3 specifies the number of 
circumferential and radial stations used in the computational mesh. Name- 
list LlST4 specifies the thermodynamic model and the molecular transport 
properties. Namelist LIST5 contains the parameters -which specify the center- 
body and the cowl contours. Namelist LIST6 contains all convergence tolerances 
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and iteration loop limits. Namelist LIST7 contains the parameters that specify 
whether or not any debug output is to be printed. 

In addition to the namelist inputs, the two-dimensional initial -value 
plane is read in by a formatted read of file ITAPE. Alternatively, if the 
forebody is conical ahead of the axial location of the initial -value plane, 
the initial -value plane data may be generated internally in the computer 
program by an approximate technique. 

All input data are read in within routine LINK! 0 . After the data have 
been entered, selected input parameters are checked for errors and consistency. 

If .any errors are detected, appropriate messages are printed and the execution 
of the program is aborted. The input parameters are tested in subroutine INITIL. 

It should be noted that many input parameters have default values and do 
not have to be specified. The default values for most of the input parameters 
are specified in LINKOO. 

After all input parameters and data have been entered, parameter initial- 
ization (for most of the subroutines) is performed in subroutine INITIL. 


5. INITIAL-VALUE PLANE 

The initial data must be specified on a space-like plane of constant x 
(see Figure 1). The flow must be supersonic at every point on this plane. 

The x-coordinate axis constitutes the longitudinal axis of the centerbody and 
the cowl. Moreover, the mean flow direction is assumed to be in the x- 
coordinate direction. 

If the forebody flow field is to be determined, the initial -value plane 
must be specified at an axial (x) station that is upstream of the forebody 
computational flow regime (see Figure 1). The solution is then found along 
the streamlines that pass through the data points specified on the initial - 
value plane, although some streamline addition and deletion are performed on 
the ensuing solution planes. The axial location of the last solution plane 
of the forebody flow field is adjusted to lie at the axial station of the cowl 
lip, and constitutes the initial -value plane for the internal flow field 
computation. 

If only the internal flow field is to’ be determined, the initial -value 
plane must be specified at the axial station which corresponds to the x-position 
of the cowl lip (see Figure 1). The cowl lip is assumed to be contained in 
a plane of constant x. For the internal flow field integration, a point re- 
distribution is performed on the initial -value plane. This point redistribution 
is required to define streamlines which lie in the stream surface formed by 
the cowl boundary. The solution is then found along the streamlines that 
pass through the redistributed points on the plane at the cowl lip axial 
station. It should be noted that the internal flow field cannot be computed 
if the bow shock wave is ingested into the annulus- 
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The initial -value plane may be specified externally by the user. The user 
supplied data is entered by a formatted read of the file ITAPE. 

The solution obtained by Jones (3) for the flow about a circular cone at 
nonzero incidence has been well substantiated. Using a conversion algorithm, 
the results of that computer program can be made compatible with the input 
data required by the computer program developed in the present investigation. 
For situations in which the forebody is conical up to the axial station where 
the initial-value plane is located, the Jones program is the recommended source 
for the initial data. 

Alternatively, if the forebody is conical ahead of the axial location 
of the initial -value plane, the flow property field on the initial -value 
plane can be generated internally in the computer program. The internally 
generated initial -value plane is obtained by an approximate technique which 
employs the Taylor-Maccoll solution for the flow about a circular cone at 
zero incidence. A superposition procedure is used to obtain an approximation 
to the flow about a circular cone at nonzero angle of attack by neglecting the 
cross flow effects. This superposition procedure effectively amounts to com- 
puting the flow turning angle in the meridional plane of the given solution 
point, and then obtaining' the flow properties at that point by applying the 
Taylor-Maccoll solution for a cone half-angle equal to the flow turning angle. 
It must be emphasized that this is only an approximate technique, giving the 
well accepted Taylor-Maccoll solution at zero incidence, but becoming in- 
creasingly less accurate as the angle of attack is increased. Subroutines 
IVSGEN , TMC0NE, INTGRT, L0GIC, C0NVRG, 0DESLV, and FUNCTM are used to generate 
internally the initial-value plane flow property field. 


6. THERMODYNAMIC MODEL 

The assumed thermodynamic model is that of a thermally and calorically 
perfect gas. Other thermodynamic models can be incorporated into the computer 
program by suitably modifying or replacing subroutines THERM0, STAG, TGRAD, 
JUMP, TMC0NE , L0GIC, and FUNCTM. 


7. MOLECULAR TRANSPORT PROPERTIES 

Both the dynamic viscosity and the thermal conductivity are assumed to be 
functions of the temperature only. The assumed functional form for the dynamic 
viscosity is given by Sutherland’s law (2). The assumed functional form for 
the thermal conductivity is a quadratic function of temperature. Other repre- 
sentations for the dynamic viscosity and the thermal conductivity can be 
incorporated into the analysis by suitably modifying or replacing subroutine 
TPR0PS . 
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8. CENTERBODY AND COWL .CONTOURS 

The computer program developed in the present investigation assumes that 
both the centerbody contour and the cowl contour are axi symmetric. For the 
purposes of geometry description, the axial (x) domain is divided into a num- 
ber of intervals. In any interval, the body radius may be specified by either 
tabular input, or by supplying the coefficients of a cubic polynomial written 
as a function of x. It is assumed that both the centerbody contour and the 
cowl contour are smooth and have continuous first derivatives. More arbitrary 
geometries, such as those having elliptic or super-elliptic cross-sections, 
can be incorporated into the analysis by suitably modifying or replacing 
subroutine BJ3UNDY. 


9. FLOW SYMMETRY 

Four flow symmetry options are incorporated into the computer program. 

The most general case is when no planes of symmetry exist. This -option is 
used to compute the flow field for fully three-dimensional inlets at in- 
cidence. The second case is when one plane of flow symmetry exists. This 
option is used for computing the flow field for axi symmetric inlets at angles 
of attack. This second case of flow symmetry is the one most likely to arise 
in the class .of problems being considered in the present investigation. The 
third case is when two planes of flow symmetry exist. This option is used to 
compute the flow field for three-dimensional inlets with two planes of geometric 
symmetry at zero angle of attack. The final option is when the flow is 
axisymmetric. This option is -used to compute the flow field in axisymmetric 
inlets at zero incidence. 

When determining the interpolation polynomial fit point stencils, solu^- 
tion point reflections must be performed when flow symmetry exists. Subroutines 
REFLEK and INTRFL are used for this purpose. 


10. OUTPUT 

Preliminary information is printed by LINK10. The initial -value plane 
and all solution planes are printed by subroutine PRN0UT. The solution points 
on the space curves defined by the intersection of the internal shock wave with 
the solid boundaries are output by subroutine SHKRFL. In addition to the 
position of and the dependent variables at a solution point, the Mach number, 
static temperature, velocity magnitude, stagnation pressure, and stagnation 
temperature are also printed. The mass flow rate across every solution plane, 
calculated by trapezoidal rule integration in subroutine MASS, is also printed. 

The quantity of results which is printed is controlled by the input para- 
meter KPRINT. Three printout options exist: body streamline points and shock 
wave points only, all solution points, and all solution points and the normal 
vector components and incident normal Mach numbers for the shock wave surface. 
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SECTION III 


SUBROUTINE DESCRIPTIONS 


1. INTRODUCTION 

» 

In this section, a brief description is given of the function of each sub- 
routine in the computer program. This information supplements the information 
available in the the form of comment statements within the program. 


2. OVERLAY (0,0} 

LINKOO . This program routine is the main control routine in 0VERLAY (0,0), 
the resident overlay. LINKOO first calls 0VERLAY {1,0} for data input, para- 
meter initialization, and, if desired, internal generation of the flow property 
field on the initial-value plane. LINKOO then calls 0VERLAY (2,0) to perform 
the flow field integration. Most of the program constants and input parameters 
have their default values specified in LINKOO. Moreover, the reader and 
printer call numbers., denoted by IRE'and IWR, respectively, are intialized 
in LINKOO. 

B0UNDY . This subroutine is used for the specification of both the fore- 
body/centerbody and cowl geometries. The version of B0UNDY supplied with the 
program assumes that both the forebody/centerbody and the cowl are axisymmetric. 
More arbitrary geometries can be described by suitably modifying or replacing 
this subroutine. If B0UNDY is replaced, the subroutine argument list must be 
identical to the existing list. The parameters in the argument list are de- 
fined as follows. 

XABS x-coordinate 

Y y-coordinate 

Z z-coordinate 

ALPHA polar angle defined by tan (z/y) 

RB0DY radius of either forebody/centerbody or cowl 

BNX x-component of outward body normal unit vector to forebody/ 
centerbody or cowl 

BNY y-component of outward body normal unit vector to forebody/ 
centerbody or cowl 

BNZ z-component of outward body normal unit vector to forebody/ 
centerbody or cowl 

J ' If J=1 ,' forebody/centerbody geometry is to be specified. 

If J=2, cowl geometry is to be specified. 
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K If K=0, read in XABS and ALPHA, and compute RB0DY and. the cor- 

responding Y and Z. If K=1 , read in XABS and ALPHA, and compute 
RB0DY, Y, Z, BNX , BNY, and BNZ. If K=2, read in XABS, Y, and Z 
(coordinates of a point not on the body), and compute RB0DY, 

ALPHA, Y, Z, BNX, BNY, and BNZ at the point on the body where 
the projection of the body normal in the (y,z)-plane passes 
through the originally specified point. 

THERMO . This subroutine computes the temperature, sonic speed., and the 
thermodynamic parameter £. The assumed thermodynamic model is that of a 
thermally and calorically perfect gas. Other thermodynamic models . may be 
incorporated into the analysis by suitably modifying or replacing this sub- 
routine. If subroutine THERM0 is replaced, the argument list must be identical 
to the existing list. The parameters in the argument list are defined as 
follows. 


P 

R0 

T 

A 

C0EFF 

KT 

KA. 

KC 


pressure (p) 
density (p) 
temperature (t) 
sonic speed (a) 
£=(l/pt)(3p/3s) p , where s 
If KT = 1, compute T. If 
If KA = 1 , compute A. If 
If KC = 1 , compute C0EFF. 


is the entropy per unit mass 
KT f 1 , do not compute T. 

KA f 1 , do not compute A. 

If KC f 1 , do not compute C0EFF. 


Thermodynamic property information is supplied to subroutine THERM0 through 
the named common block GAS! . This common block conta-ins the following para- 
meters: RUNIV (universal gas constant), TDT (temperature data array), WTM0L 
(molecular weight data array) , CP (specific heat data array), NGAS1 (an integer 
to denote the number of elements in the TDT, WTM0L, and CP arrays), and KGAS’l 
(an integer to denote which gas model is to be used). These parameters are 
included for modifying subroutine THERM0 to include real gas effects. 

STAG. This subroutine computes the Mach number, the stagnation pressure, 
and the stagnation temperature. The assumed thermodynamic model is that of a, 
thermally and calorically perfect gas. Other thermodynamic models may he in- 
corporated into the analysis by suitably modifying or replacing this subroutine. 
If STAG is replaced, the argument list must be identical to the existing list. 
The parameters in the argument list are defined as follows. 

P static pressure 

R0 static density 

T static temperature 

A sonic speed 

0 velocity magnitude 

XM Mach number 
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PT stagnation pressure 

TT stagnation temperature; 

Thermodynamic property information is supplied to subroutine STAG through 
the named common block GAS! . 

TGRAD . This subroutine computes the spatial gradients of temperature 
by using an analytically differentiated form of the thermal equation of state 
so that temperature derivatives can be expressed in terms of the derivatives 
of pressure and density.. The assumed thermodynamic model is that of a ther- 
mally perfect gas. Other thermodynamic models may be incorporated into the 
analysis by suitably modifying or replacing this subroutine. If TGRAD is 
replaced, the argument list must be identical to the existing list. The 
parameters in the argument list are defined as follows. 

T temperature 

P pressure 

R0 density 

DPX,DPYjDPZ first partial derivatives of pressure with respect to 

x, y, and z, respectively 

DPXX>DPYY,DPZZ second partial derivatives of pressure with respect 

to x, y, and z, respectively 

DR0X,DR0Y,DR0Z first partial derivatives of density with respect 

to x, y, and z, respectively 

DR0XX,DR0YY,DR0ZZ second partial derivatives of density with respect 

to x, y, and z, respectively 

DTX,DTY,DTZ first partial derivatives of temperature with 

respect to x, y, and z, respectively 

DTXX,DTYY ,DTZZ second partial derivatives of temperature with 

respect to x, y, and z, respectively 

Thermodynamic property information is supplied to subroutine TGRAD 
through the named common block GAS1 . 

TPR0PS . This subroutine computes the dynamic viscosity, the thermal 

conductivity, and their spatial gradients. The assumed functional form for 
dynamic viscosity is given by Sutherland's Law (2). The assumed functional 
form for thermal conductivity is a quadratic curve fit as a function of tem- 
perature. Other transport property formulations may be incorporated into the 
analysis by suitably modifying or replacing this subroutine. If TPR0PS is 
replaced, the argument list must be identical to the existing list. The 
parameters in the argument list are defined as follows. 

T static temperature 

p static pressure 

VIS dynamic viscosity 
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C0N thermal conductivity 

DTX»DTY>DTZ first partial derivatives of temperature with respect 

to x s y, and 2 , respectively 


DPX»DPY>DPZ first partial derivatives of pressure with respect to 

x, y, and z, respectively 

DVISX,DVI$Y,DVISZ first partial derivatives of viscosity with respect 

to x, y, and z, respectively 

DC0NX,DC0NY,DC0NZ first partial derivatives of thermal conductivity 

with respect to x, y, and z, respectively 


Transport property information is ‘supplied to subroutine TPR0PS through 
the named common block GAS2. This common block contains the following para- 
meters: TDL (temperature data array), VISD (viscosity data array), C0ND 
(thermal conductivity data array) , NGAS2 (an integer to denote the number 
of elements in the TDL, VISD, and C0ND arrays), and .KGAS2 (an integer to 
denote which transport property model is to be used). These parameters are 
included for modifying TPR0PS to use tabular data. 

JUMP . This subroutine employs the Rankine-Hugoniot relations to compute 
the downstream flow properties at a shock wave solution point. The Shock wave 
surface normal unit vector components and the flow properties upstream of the 
shock wave are entered in the argument list. The computed downstream flow 
properties are also transmitted in the argument list. The assumed thermodynami 
model is that of a thermally and calorically perfect gas. Other thermodynamic 
models may be incorporated into the analysis by suitably modifying this sub^- 
routine. Thermodynamic property information is supplied to subroutine JUMP 
through the named common block GAS1. 


GELG. This subroutine is used for solving a system of simultaneous linear 
equations. The system is solved using Gaussian elimination with complete 
pivoting. This subroutine is IBM library subroutine GELG. 


3. 0VERLAY (1,0) 

LINK10 . This program routine is the main control routine in 0VERLAY (1,0) 
All input parameters and initial data are entered in LINK! 0 . In all cases, 
the seven namelists LIST! to LIST7 are entered. If the flow property field 
on the initial -value plane is externally specified, a formatted read of file 
ITAPE is used to enter this information. After the input data have been 
entered, subroutine INITIL is called for testing for input errors and per- 
forming parameter initialization. If the flow property field on the initial- 
value plane is to be internally generated, subroutine IVSGEN is then called 
from subroutine INITIL. After subroutine INITIL has been called, preliminary 
output is printed by LINK10. 
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INITIL . This subroutine is called from routine LINK10, and is used to 
test for errors on selected input parameters and to perform parameter initial- 
ization. If an input parameter error is found, an appropriate message is 
generated and the program execution is aborted. After the selected input 
parameters have been tested, parameter initialization is performed. If the 
flow property field on the initial -value plane is to be internally generated, 
subroutine IVSGEN is then called. 

IVSGEN. This subroutine is used in conjunction with subroutines TMC0NE, 
INTGRT, L0GIC, C0NVRG , 0DESLV, and FUNCTM to internally generate the flow 
property field on the initial -value plane if- the forebody is conical ahead of 
the initial -value plane. The internally generated initial -value plane is ob- 
tained by an approximate technique which employs the Taylor-Maccoll solution 
for the flow about a circular cone at zero incidence. A superposition pro- 
cedure is used to obtain an approximation to .the flow about a circular cone at 
nonzero angle of attack by neglecting the cross flow effects. This super- 
position procedure effectively amounts to computing the flow turning angle in 
the meridional plane of the given solution point, and then obtaining the flow 
properties at that point by applying the Taylor-Maccoll solution for a cone 
half-angle equal to the flow turning angle. The shock wave angle is then 
measured from the original streamline direction in the appropriate meridional 
plane. It must be emphasized that this is only an approximate technique, 
giving the well accepted Taylor-Maccoll solution at zero angle of attack, but 
becoming increasingly less accurate as the angle of attack is increased. 

TMC0NE . This subroutine is one of the subroutines used for the internal 
generation of the initial -value plane. TMC0NE is called from subroutine 
IVSGEN, and performs initialization and post calculation storage of the 
Taylor-Maccoll solution used in computing the initial -value plane flow pro- 
perty field. The assumed thermodynamic model is that of a thermally and 
calorically perfect gas. Other thermodynamic models may be incorporated into 
the analysis by suitably modifying this subroutine. Thermodynamic property 
information is supplied to subroutine TMC0NE through the named corrcnon block 
GAS1. 


INTGRT . This subroutine is one of the subroutines used in internally 
generating the flow property field on the initial -value plane. INTGRT in- 
tegrates the Taylor-Maccoll equation inward towards the forebody until a point 
is reached where the normal component of velocity is sufficiently small. 

L0GIC . This subroutine is one of the subroutines used in internally 
generating the flow property field on the initial -value plane. L0GIC is used 
to control- the computation of the Taylor-Maccoll solution in a particular 
meridional plane. The assumed thermodynamic model is that of a thermally and 
calorically perfect gas. Other thermodynamic models may be incorporated into 
the analysis by suitably modifying this subroutine. Thermodynamic property 
information is supplied to subroutine L0GIC through the named common block 
GAS! . 
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C0NVRG . This subroutine is one of the subroutines used in internally 
generating the flow property field on the initial -value plane. C0NVRG is 
used in testing if the normal component of velocity has vanished at the body 
in the external flow field integration. 

0DESLV. This subroutine is one of the subroutines used in internally 
generating -the flow property field on the initial -value plane. 0DESLV in- 
tegrates the Taylor-Maccol 1 equation by employing a fourth-order Ru'nge-Kutta 
method. 

FUNCTM . This function is used in internally generating the flow property 
field on the initial-value plane. FUNCTM evaluates the Taylor-Maccol 1 equation. 
The assumed thermodynamic model is that of a thermally and calorically perfect 
gas. Other thermodynamic models may be incorporated into the analysis by 
suitably modifying this function. Thermodynamic property information is 
supplied to FUNCTM through the named common block GAS1. 


4. 0VERLAY (2,0) 

LINK20 . This program routine is the main control routine in 0VERLAY (2,0). 
LINK20 calls the integration control overlays, 0VERLAY (2,1), 0VERLAY (2,2), 
and 0VERLAY (2,3). If a program restart is specified, the last secondary level 
overlay that was in use at the time the initial execution was terminated is 
called. 

REDIST . This subroutine is used to perform point redistribution on the 
solution plane at the axial station of the cowl lip. This point redistribution 
is required to obtain a uniform point distribution and to obtain streamlines 
tfhich lie in the stream surface formed by the cowl boundary. The redistributed 
points are arranged symmetrically in the computed sector. These points lie 
pn rays which have equal angular increments from one another, with the points 
pn each ray being spaced at equal radial increments. The flow properties 
it these points are obtained by bivariate interpolation. If the viscous flow 
pption is specified and the internal flow field integration option in which 
shock waves are not discretely fitted is specified, point redistribution is 
ilso performed on the forebody flow field solution plane immediately upstream 
pf the solution plane located at the cowl lip axial station. 

STNSEL . This subroutine is employed to determine the solution point on a 
jiven solution plane that is nearest to an arbitrary point whose coordinates 
ire supplied in the argument list. After the proper solution point has been 
ietermined by a search of the computed sector, it is used as the base point 
)f the nine point stencil used in formulating the quadratic bivariate inter- 
iolation polynomial. STNSEL calls subroutine BIPTS in the course of determining 
:he proper solution point and in formulating the interpolation polynomial 
:oefficients. 
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REFLEK . This subroutine performs reflections, about the coordinate axes, 
of the solution points in the computed sector on a given solution plane when 
flow symmetry exists. The solution point reflections are required to obtain 
the appropriate fit point stencils used in formulating the univariate, bi- 
variate, and trivariate interpolation polynomials. REFLEK is used to reflect 
points for the case of one plane of flow symmetry, the case of two planes of 
flow symmetry, and the axisymmetric flow case. Both streamline and shock wave 
points are reflected in REFLEK. 

XSTEP . This subroutine determines the axial marching step allowed by the 
Courant-Friedrichs-Lewy (CFL) stability criterion. Except in the vicinity of 
an internal shock wave intersection with a solid boundary, the marching step 
computed by XSTEP is that used to regulate the distance between successive 
solution planes. In the vicinity of an internal shock wave-solid boundary 
intersection, special constraints are used to determine the marching step. 

XSTEP applies the CFL stability criterion to all streamline points in the 
computed sector. The actual marching step is taken as the allowable x-step 
at the most restrictive point. The stability criterion i,s not applied to the 
shock wave points. Moreover, the internal flow field shock wave points are 
ignored in defining the convex hull of the finite difference network when 
application of the stability criterion is made to a streamline point. 

MASS . This subroutine is used to compute the mass flow rate across a given 
solution plane. The mass flow rate is numerically ascertained by a trapezoidal 
rule integration procedure in which the incremental mass flow rate is computed 
through an elemental triangle formed by three adjacent points. The total of 
these incremental mass flow rates is the total mass flow rate. If flow sym- 
metry exists, only the mass flow rate in the computed sector is found by 
integration. The mass flow rate across the entire plane is then obtained by 
use of an appropriate multiplier. Special logic is used to trace the internal 
shock wave in the integration procedure when the mass flow rate across an 
internal flow field solution plane is being computed. 

PRN0UT . This subroutine is used to print the integration results for 
all solution planes. Both the external flow field and the internal flow 
field solution planes are printed by PRN0UT. In addition, the initial-value 
plane, the redistributed data plane at the cowl lip axial station, and the 
restart plane are output by PRN0UT. Three print options (specified .in the in- 
put by the input parameter KPRINT) are available: body streamline points and 
shock wave points, all solution points, and all solution points and shock 
wave parameters (incident normal Mach number and shock wave surface normal 
unit vector components). The solution points along the space curves formed by 
the intersections of the internal shock wave with the solid boundaries are not 
printed by subroutine PRN0UT. They are printed by subroutine SHKRFL. 

UNIPTS . This subroutine selects the fit points used in formulating the 
univariate interpolation polynomials that are used to describe the shock wave 
radius and the shock wave angle along the curve defined by the intersection 
of the shock wave with a given solution plane. Three adjacent shock wave 
solution points constitute the fit point stencil in regions that are away 
from an internal shock wave-solid boundary intersection. In the region of 
an intersection, the fit point stencil may be appropriately expanded to be 
in accord with the Courant-Friedrichs-Lewy (CFL) stability criterion. 
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UNI FIT . This subroutine is used to determine the coefficients in the 
quadratic univariate interpolation polynomial. The three coefficients in this 
polynomial are determined by a fit to three data points which are supplied 
to UNIFIT through the call statement. Subroutine GELG is called to solve the 
system of simultaneous linear equations that determines the coefficients. 

BIPTS . This subroutine selects the fit points used in formulating the 
quadratic bivariate interpolation polynomials that are used for intraplanar 
flow property determination. A base point and its eight immediate neighbors 
constitute the fit point array. Two types of fit point stencils are used: 
interior point and boundary point. A boundary point stencil is employed when 
the interpolation base point (the fit point nearest to the interpolated point) 
is on the shock wave. Special logic is used to insure that no point stencil 
bridges the shock wave. 

BIF1T . This subroutine is used to determine the coefficients in the 
quadratic bivariate interpolation polynomial. The six coefficients in this , 
polynomial are determined by a least squares fit of nine data points which are 
supplied to BIFIT through the named common block FITPTS. Subroutine SYMSIM 
is called to solve the system of simultaneous linear equations (normal equa- 
tions) which determines the coefficients in the interpolation polynomial. 

This system of linear equations has a symmetric coefficient matrix. 

TRIPTS . This subroutine selects the fit points used in formulating both 
the linear trivariate interpolation polynomial and the quadratic trivariate 
interpolation polynomial. Four solution points are used in formulating the 
linear trivariate interpolation polynomial. Fourteen solution points are used 
in formulating the quadratic trivariate interpolation polynomial. The linear 
polynomial is used for flow property interpolation on the upstream side of the, 
shock wave surface. The quadratic polynomial is used for flow property inter- 
polation on both the downstream side of the shock wave surface and on the 
stream surface formed by a solid boundary. 

TRIFTU . This subroutine is used to determine the coefficients in the 
linear trivariate interpolation polynomial. The four coefficients in this 
polynomial are determined by a fit to four data points which are supplied to 
TRIFTU through the named common block FITPTS. Subroutine GELG is called to 
solve the system of simultaneous linear equations which determine the coef- 
ficients. 

TRIFTD . This subroutine is used to determine the coefficients in the 
quadratic trivariate interpolation polynomial. The eight coefficients in this 
polynomial are determined by a least squares fit of fourteen data points which 
are supplied to TRIFTD through the named common block FITPTS. Subroutine 
SYMSIM is called to solve the system of simultaneous linear equations (norma] 
equations) which determines the coefficients. This system of equations has a 
symmetric coefficient matrix. 

SYMSIM . This subroutine solves a system of simultaneous linear equations 
with a symmetric coefficient matrix. The system is solved using Gaussian 
elimination with pivoting in the main diagonal. This subroutine is a modified 
version of IBM library subroutine GELS. 
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SHKINT. This subroutine is used to evaluate the univariate interpolation 
polynomials that are employed in describing the shock wave radius and the 
shock wave angle along the curve that is defined. by the intersection of the 
shock wave with a given solution plane. The interpolation polynomial coef- 
ficients are supplied to SHKINT through the named common block SINTRP. The 
independent variable is the polar angle THETA, which is transmitted through 
the subroutine call statement. 

PRPINT . This subroutine is used to evaluate the bivariate and tri variate 
interpolation polynomials that are used for flow property determination on 
solution planes, shock waves, and solid boundary stream surfaces. The inter- 
polation polynomial coefficients are supplied to PRPINT through the named 
common block INTRP. The independent variables are the three position coordi- 
nates XA, Y, and Z, which are transmitted through the subroutine call 
statement. 

FDERIV . This subroutine is used to compute the first partial derivatives 
of the flow properties on a solution plane, a shock wave, or a solid boundary 
stream surface. The derivatives are obtained by evaluating an analytically 
differentiated form of the appropriate interpolation polynomial. The inter- 
polation polynomial coefficients are transmitted to FDERIV through the named 
common block INTRP. The independent variables are the three position coordi- ■ 
nates XA, Y, and Z, which are transmitted through the subroutine call 
statement . 

VECTOR . This subroutine is used in conjunction with subroutine SH0CK 
to determine the components of the unit vector 0 which is used in the para- 
meterization of the wave surface compatibility relation. The unit vector 0 
is orthogonal to the velocity vector that is downstream of the shock wave at 
the shock wave solution point, and has its projection on the (y,z)-plane lie 
in the meridional plane which passes through the shock wave solution point. 

PLANAR . This subroutine is used to compute the parameters employed in 
the formulation that represents the shock wave surface. This formulation is 
then used in determining the intersection point of either a streamline or a 
bicharacteristic with the shock wave surface. Additionally, PLANAR initializes 
some parameters used in determining the intersection point of a bi characteristic 
with the stream surface formed by a solid boundary. 

INTSCT . This subroutine is used to compute the intersection point of 
either a streamline with the' shock wave surface, or the intersection point of 
a bicharacteristic with either the shock wave surface or the stream surface 
formed by a solid boundary. The determination of the intersection point 
coordinates is performed in an iteration loop which uses the secant method to 
relax the difference in the radius of the point of intersection obtained from 
integration of the equation for a streamline or bi characteristic and that 
obtained from the appropriate surface formulation. 

FUNCT . This subroutine is used in conjunction with subroutine INTSCT to 
determine a streamline-shock wave intersection point, or the intersection point 
of a bicharacteristic with either a shock wave or a solid boundary. The axial 
position of the assumed intersection point is supplied to FUNCT through the 
subroutine call statement. From the given axial position, FUNCT computes the 
corresponding y and z coordinates of the assumed intersection point by both 
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integrating the equation for a streamline or bicharacteristic, and by eval- 
uating the appropriate surface formulation. The difference in the radius 
obtained from the streamline or bicharacteristic equation integration and that 
obtained from the surface formulation is reduced to within a specified 
tolerance of zero by the iteration, method used in subroutine INTSCT. 

RADIUS . This subroutine is used to compute the radius of a point on the 
shock wave surface or on the stream surface formed by a solid boundary. The 
axial position and the polar angle of the point are supplied to RADIUS through 
the subroutine call statement. To obtain the shock wave radius at the- desired 
point, a linear interpolation is performed in the meridional plane of the 
point between two space curves which are defined by either a shock-wave solution 
plane 'intersection or by a shock wave-solid boundary intersection. To obtain 
the body radius at the desired point, subroutine B0UNDY is called. 

XFIT . This subroutine is used to curve fit, as a function of the polar 
angle, the axial (x) position of an internal shock wave-solid boundary inter- 
section. A quadratic polynomial expressed in terms of the polar angle is used 
for this representation. 

FRCFNS . This subroutine is used to compute the molecular transport 
forcing terms used in both the governing equations and the compatibility re- 
lations. FRCFNS is called by subroutines S0LVE and SH0CK. The molecular 
transport terms can be included in the computation of the external flow field 
about the forebody or in the computation of the internal flow field in which 
shock waves are not discretely fitted. The program option in which internal 
shock waves are discretely fitted does not have the capability to include the 
influence of molecular diffusion in the computation, but rather assumes the 
flow to be inviscid and adiabatic. 

S0LN . This subroutine calls subroutines L0CATE and S0LVE for computing 
either a solid boundary solution point or an interior solution point. 

L0CATE . This subroutine is used to compute the locations of and the 
flow properties at the streamline and bicharacteristic intersection points for 
the standard interior point, standard solid boundary point, shock-modified 
interior point, and shock-modified solid boundary point unit processes. The 
point locations and flow properties at these points are transmitted to sub- 
routine S0LVE through the named common blocks RELAY! and RELAY2. Subroutine 
S0LVE then solves the system of compatibility relations to obtain the position 
of and the flow properties at the solution point. 

S0LVE . This subroutine is used to solve the system of compatibility 
equations for the standard interior point, standard solid boundary point, 
shock-modified interior point, and shock-modified solid boundary point unit 
processes. If the viscous and thermal diffusion terms are to be included in 
the computation, S0LVE calls subroutine FRCFNS. The system of five compat- 
ibility relations is solved by calling subroutine GELG. ’ 
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SH0CK . This subroutine is used to compute the solution for all field- 
shock wave points. For the bow shock wave points, the free-stream flow 
conditions are used for the upstream flow properties. For the internal flow 
field shock wave points, an interior point or solid boundary unit process is 
applied to obtain the upstream flow properties. Body solution points on the 
downstream side of the cowl lip shock wave or on the downstream side of a 
reflected internal shock wave are computed using the solid boundary- shock wave 
point unit process (subroutine BSH0CK) . 


5. 0VERLAY (2,1) 

LINK21 . This program routine contains the control logic used in the 
computation of the external flow field about the forebody. LINK21 is the 
main control routine in 0VERLAY (2,1). 0VERLAY (2,1) is a secondary level 
overlay which is called from routine LINK20. 

RCNTRL . This subroutine is used to control the number of radial stations 
on successive solution planes in the forebody flow field integration. If 
a sufficient influx of mass across the bow' shock wave has occured between the 
current solution plane and the last solution plane where point addition or 
deletion was performed, a new ring of interior field points is added between 
the ring of shock wave points and the outermost ring of existing interior 
field points. If, after a number of successive point additions, a specified 
number of radial stations has been reached, point deletion is performed. Here, 
selected interior field points are deleted from the storage arrays while the 
bow shock wave points and the body streamline points are retained. 


6. 0VERLAY (2,2) 

LINK22 . This program routine contains the control logic used in the com- 
putation of the internal flow field in which shock waves are discretely fitted. 
LINK22 is the main control routine in 0VERLAY (2,2). 0VERLAY (2,2) is a 
secondary level overlay which is called from routine LINK20. 

BSH0CK. This subroutine is used to compute the flow properties ata point 
on the body that is downstream of either the cowl lip shock wave or an internal 
reflected shock wave (solid body-shock wave point unit process). This sub- 
routine is used only in the option where the internal flow field integration 
in which shock waves are discretely fitted is employed. 

SPRJCT . This subroutine is used in the internal flow field integration 
to project the internal shock wave from the current initial -value plane to 
the current solution plane. The projected intersection of the internal shock 
wave with the solution plane is then used to determine which streamlines do 
and do not penetrate the shock wave. 

PENTRE . This subroutine is used in the internal flow field integration 
to control the computation of an interior field point when the streamline has 
penetrated the internal shock wave. If the streamline-shock wave intersection 
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point is sufficiently close to the current solution plane, an interior point 
unit process on the downstream side of the shock wave- is not performed. In- 
stead, in this case, a streamline projection onto the solution plane and 
subsequent flow property interpolation in this plane is performed. Alterna- 
tively, if the streamline-shock wave intersection point is sufficiently far 
from the solution plane, an interior point unit process is performed on the 
downstream side of the shock wave. 

TRACE . This subroutine is used in conduction with subroutine INTSCT to 
compute the intersection point of a streamline with an internal shock wave. 

With the intersection point coordinates determined, the trivariate interpolation 
polynomials are evaluated to obtain the flow properties at the- intersection 
point. 

S0LINT . This subroutine is used in the internal flow field integration 
to determine the intersection point coordinates of and the flow properties at 
an interior field point when the streamline penetrates an internal shock wave 
with the intersection point being sufficiently close to the current solution 
plane. S0LINT is called from subroutine PENTRE. 

SWITCH ., This subroutine is used in the internal flow field integration 
to perform the post computation interchange of indices when a streamline which 
initially appeared to intersect the internal shock wave ultimately did not. 
Consequently, a standard interior point unit process is used for this point 
instead of a shock-modified interior point unit process performed on the 
downstream side of the shock wave. 

SHKRFL . This subroutine contains the control logic used for calculating 
an internal shock wave-solid boundary intersection. 

INTRFL , This subroutine performs the necessary point reflections when 
flow symmetry exists for the solution points on the space curve defined by 
the intersection of the internal shock wave with a solid boundary. 

ISH0CK . This subroutine is used in the internal flow field integration 
to compute the incident shock wave upstream and downstream flow properties 
at a point where the incident shock wave intersects a solid boundary. 

STRSHK . This subroutine is used to compute the position of and the flow 
properties at the intersection point of a body streamline with the space curve 
defined by the intersection of the incident shock wave with a solid boundary. 

RSH0CK. This subroutine is used with subroutine BSH0CK to compute the 
flow properties on the body downstream of an internal reflected shock wave. 


7. 0VERLAY (2,3) 

LINK23. This program routine contains the control logic used in the com- 
putation of the internal flow field in which shock waves are not discretely 
fitted. LINK23 is the main control routine in 0VERLAY (2,3). 0VERLAY (2,3) 
is a secondary level overlay which is called from routine LINK20. 
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SECTION IV 


INPUT PARAMETERS 


1 . INTRODUCTION 

The input data required for execution of the computer program are enteVed 
by both namelist input and formatted r.ead statements. In all cases, the seven 
namelists LIST1 , LIST2 , LIST3 , LIST4, LIST5, LIST6 , and LIST7 are entered.* In 
cases in which the user selects to specify the flow property field on the 
initial-value plane, that information is entered by a formatted read of either 
TAPE 5 (the input file) or TAPE 3. 

In general, only those parameters and data pertinent to the particular 
problem being considered must be entered. Many input parameters have default 
values and do not need to be specified unless values other than the default 
values are to be entered. 

In this section, each input parameter is defined. Where applicable, both 
the default value and the typical value of the -input parameter are given. 


2. TITLE CARD 

The first card of each data deck is a title card on which 72 'al phanumeric 
characters (any standard Fortran characters) of identifying information may 
be specified. This card must be the first card of the data deck even if no 
information is listed on it. The format of the- card is (12A6). 


3. NAMELIST LIST! 

.The parameters entered in namelist LIST1 control the overall execution of 
the program. 

KUNIT An integer variable denoting whether English absolute units or SI 
units are to be used in the computation. If KUNIT = 1, English 
absolute- units are employed. If KUN-IT = 0, SI units are employed. 

A default value of 1 is specified for KUNIT. 

KCALL A one-dimensional integer variable array consisting of three ele- 
ments. Each element of KCALL specifies whether or not a particular 
flow field integration option is to be performed. If KCALL(I) = 1- 
(1=1 ,2,3), then the corresponding flow field integration option is 
performed. If KCALL(I) = 0, the integration option is not performed 
The elements of KCALL control the flow field integration options and 
have specified default- values as listed, below. 
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Default Value 


KCALL(I) 



KCALL(l) forebody flow field 

KCALL (2 } internal flow field with 

discrete shock wave fitting 

KCALL (3) internal flow field without 

discrete shock wave fitting 


1 

1 

0 


Specifying KCALL(2) = 1 and KCALL(3) = 1 simultaneously wi 11 cause 
an error message to be printed and the program execution to be 
aborted. 


XEND A one-dimensional real variable array consisting of three elements. 

XEND(I) (1=1 ,2 ,3) denotes the x-position, in either feet or meters, 
at which the flow field integration specified by the corresponding 
element of KCALL is to be terminated. Each element of XEND denotes 
a flow field integration option termination point and has a default 
value as follows. 


XEND ( I ) 

Termination Point for 

Default Value 

XEND(l) 

forebody flow field 
itegration 

2.0 ft 

XEND(2) 

internal flow field 
integration with discrete 
shock wave fitting 

3.5 ft 

XEND ( 3 ) 

internal flow field 

3.5 ft 


integration without discrete 
shock wave fitting 



If any element of XEND exceeds the x-position to which the center- 
body geometry is specified, or if XEND ( 2 ) or XEND(3) exceeds the 
x-position to which the cowl geometry is specified, appropriate 
error messages are printed and the program execution is aborted. 
Each element of XEND must be positive. It should be noted that if 
KCALL(I) = 0, XEND(I) does not have to be specified. 

RCAVG A positive real variable denoting the estimated average radius, in 
either feet or meters, of the bow shock wave at XEND(l). The 
specified value of RCAVG is used in estimati ng' the captured mass 
flow rate at XEND(l). This mass flow rate is used in determining 
if point addition is to be performed in the forebody flow field 
integration. RCAVG must be specified only if the forebody flow 
field integration option .is used [KCALL(l) = 1]- RCAVG must be 
specified even if the forebody integration option is the only in- 
tegration option used [KCALL (2 ) = KCALL(3) = 0]. The cowl lip 
radius may be used for RCAVG. A default value of 0.8 ft is 
specified for RCAVG. 
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KVISCY An integer variable denoting whether or not the viscous and thermal 
diffusion terms are to be included in the computation of the fore- 
body flow field or the internal flow field in which shock waves are 
not discretely fitted. If KVISCY = 1, these terms are included in 
the computation. If KVISCY = 0, these terms are not included in 
the computation. A default value of 0 is specified for KVISCY. 

KSYM An integer variable denoting the flow symmetry option to be employed 
in the computation. KSYM can have the values 0, 1, 2 , and 3 cor- 
responding to: 


KSYM Flow Symmetry Option 

0 no planes of symmetry - computed sector is the entire 
solution plane 

1 one plane of symmetry - computed sector is the half- 
plane bounded by the y-axis and containing the +z- 
axis 

2 two planes of symmetry - computed sector is the 
quadrant bounded by the +y-axis and the +z-axis 

3 ax i symmetric flow - computed sector is the single 
circumferential station on the +y-axis 

The point networks for the four flow symmetry options are illustrated 
in Figure 3. A default value of 3 is specified for KSYM. 

KSGL0B An integer variable denoting whether or not global correction is to 
be performed in obtaining the solution for the bow shock wave points. 
If KSGL0B = 1, global correction is performed. If KSGL0B = 0, global 
correction is not performed. Global correction can only be performed 
for the bow shock wave points and not for the internal flow field 
shock wave points. Hence, KSGL0B must be specified only if 
KCALL(l) =1. A default value of 1 is specified for KSGL0B. 

KPRINT An integer variable denoting which of three print options is to be 
employed in the program execution. KPRINT can have the following 
values. 

KPRINT Print Option 

0 print body solution points and shock wave solution 
points only 

1 print all solution points 

2 print all solution points and shock wave parameters 
(incident normal Mach number and shock wave surface 
normal unit vector components) 

A default value of 1 is specified for KPRINT. 


I 
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IPRSTP A positive integer variable denoting the plane number at which the 
program execution is to be terminated. This input parameter is 
typically employed when the execution is intended to be restarted 
at plane number (IPRSTP + 1). Specifying IPRSTP £0 has no: effect 
on the execution of the program. A default value of 0 is specified 
for IPRSTP. 

KSTART An integer variable used in controlling the file operations involved 
in the restarting of the program (if desired). KSTART controls the 
storage and retrieval of the restart file as follows. 

KSTART Control of Restart File 

0 no file operations 

1 write restart file on TAPE 4 

2 read information for program restart from TAPE 4, 

and write ensuing solution planes on this tape 

A default value of 0 is specified for KSTART. TAPE 4 is linked to 
the dummy file RESTRT in the PROGRAM card. 


4. NAMELIST LIST2 

The parameters entered in namelist LIST2 specify the free-stream condi- 
tions, the inlet orientation, and the parameters which control the internal 

generation of the flow property field on the initial -val ue plane. 

MFS A positive real variable denoting the free-stream Mach number. The 

specified value of MFS must be greater than 1.0. A default value of 
3.0 is specified for MFS. 

PFS A positive real variable denoting the free-stream absolute pressure, 

in either (lbf/ft 2 ) or (N/m 2 ). A default value of 242.2 (lbf/ft 2 ) is 
specified for PFS (this value is the - pressure of the standard atmos- 
phere at an altitude of 50,000 ft). 

R0FS A positive real variable denoting the free-stream density, in either 
(slug/ft 3 ) or (kg/m 3 ). A default value of 0.0003622 (slug/ft 3 ) is 
specified for R0FS (this value is the density of the standard atmos- 
phere at an altitude ,of 50,000 ft). 

PITCH A, real variable denoting the angle, in degrees, subtended by the f.ree- 
stream velocity vector and the projection of the free-stream velocity 
vector on the (x,z)-plane, as illustrated in Figure 4. A default 
value of 0.0 degrees is specified for PITCH. 

YAW A real variable denoting the angle, in degrees, subtended by the ,x- 

axis and the projection of the free-stream velocity vector on the 
(x,z)-plane, as illustrated in -Figure 4. A default value of 0.0 
degrees is specified for YAW. 
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Vj cos (PITCH) cos (YAW) 


Vj sin (PITCH) 


Vj cos (PITCH) sin (YAW) 



XI 


A positive real variable denoting the axial (x) position, in either 
feet or meters, of the initial -value plane. If the forebody flow 
field integration option is specified [KCALL(l) = 1], X-I must be 
specified at the beginning of the forebody flow field computational 
regime (see Figure 1). If only the internal flow field integration 
option is specified [KCALL(l) = 0, KCALL(2) = 1 or KCALL(3) =1], 

XI must be specified at the axial station of the cowl lip (see 
Figure 1). XI must not fall outside of the range of axial stations 
for which the centerbody geometry is specified. Also, XI must not 
be greater than the axial station up to which the cowl geometry is . 
specified. A default value of 1.0 ft is specified for XI. 

■KIVS An integer variable denoting whether the initial -value plane is to 
be generated internally by the approximate technique described in 
Section II, or read in. If KIVS = 1, the initial-value plane is 
computed internally. If KIVS = 0, the initial-value plane must be 
supplied by the user through a formatted read of file ITAPE. The 
formatted read input is described at the end of this section. The 
internally generated initial -value plane option is applicable only 
to cases where the forebody is conical up to the axial station where 
the initial -val ue plane is located. Specifying KIVS = 1 requires 
that KSYM = 1,2, or 3and that YAW = 0.0. Furthermore, if KIVS = 1 
and KSYM = 3, PITCH must be specified as 0.0. A default value of 1 
is specified for KIVS. 

KC0N An integer variable denoting whether or not the bow shock wave is 
conical. KC0N must be specified only if the forebody flow field 
integration option is employed [ KCALL ( 1 ) = 1].. If.KC0N = 1 , the bow 
shock wave is conical. In this case, the angle at 'each initial- 
value plane bow shock wave point that is subtended by the shock wave 
and the x-axis in the meridional plane defined by the shock wave 
point is computed internally. If KC0N = 0, the bow shock wave is 
not conical, and the shock wave angles must be supplied by the user 
through a formatted read of .file ITAPE. The- formatted read input 
is described at the end of this section. If the initial-value plane 
is generated internally (KIVS = 1), a conical- bow shock wave is 
assumed (KC0N =1). A default value of 1 is specified for KC0N. 

ITAPE A positive integer variable denoting the tape number from which the 
user supplied initial -value plane is to be entered by a formatted 
■read. ITAPE must be specified only if the initial -value plane flow 
property field is to be supplied by the user (KIVS = 0). The default 

value assigned to ITAPE is 5 (the input -file). The user may specify 
ITAPE = 3, in which case the initial-value plane is read from TAPE 3. 
TAPE 3 is linked to the dummy file IVS in the PR0GRAM card. 

5. NAMELIST LIST3 

The parameters entered in namelist LIST3 specify the number of circum- 
ferential and radial stations used in the computational point network. 
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IST0P 


A positive integer variable denoting the number of circumferential 
stations used in the computed flow field sector. The value specified 
for IST0P must correspond to the flow symmetry option specified by 
KSYM (see namelist LIST1) as follows (see Figure 3). 

KSYM Allowable Value(s) of IST0P 

0 5 < IST0P < 40 

1 4 < IST0P £ 21 

2 3 < IST0P < 11 

3 IST0P = 1 

A default value of 1 is specified for IST0P [this value corresponds 
to KSYM = 3 (ax i symmetric flow)]. A value of IST0P.= 21 is recom- 
mended for KSYM = 1 . 

JMAXI A positive integer variable denoting the number of radial stations on 
the initial-value plane. Note that the initial -value plane input 
format is the same whether the computation is being started at an 
axial station upstream of the forebody computational flow regime or 
at the axial station of the cowl lip. The specified value for JMAXI 
must be at least 3 and no greater than 21. The default and recom- 
mended value for JMAXI is 11 . 

Whatever the value of JMAXI, the outermost radial station corresponds 
to the downstream bow shock wave points, and the remaining (JMAXI - 1) 
radial stations correspond to the streamline points. 

JINLET A positive integer variable denoting the number of radial stations on 
each solution plane in the internal flow field integration. JINLET 
must be specified only if the internal flow field is to be computed. 
The specified value of JINLET must be at least 4 and no greater than 
21. The recommended value and the default value of JINLET is 11. It 
should be noted that JINLET is independent of JMAXI. 

For the internal flow field integration option in which shock waves 
are not discretely fitted [KCALL(2) = 0, KCALL(3) = 1], JINLET 
specifies the number of streamline points at each circumferential 
station. For the internal flow field integration option in which 
shock waves are discretely fitted [KCALL(2) = 1, KCALL(3) = 0], the 
number of streamline points at each circumferential station is equal 
to (JINLET - 2). The remaining two storage locations are assigned 
to the upstream and downstream shock wave solution points, 

JLIMIT A one-dimensional integer variable array consisting of two elements. 

JLIMIT must be specified only if the forebody flow field integration 
option is employed [KCALL(l) = 1], The elements of JLIMIT are used 
in controlling the number of interior field points which are added 
in the computation of the forebody flow field. The first element of 
JLIMIT [JLIMIT(l)] denotes the allowable maximum number of radial 


31 



stations on a solution plane in the forebody flow field computation 
when the mass flow rate across that plane is less than a specified 
fraction [denoted by CRIT ( 7 ) in namelist LISTS] of the estimated 
mass flow rate at XEND(l), the axial location at which the forebody 
flow field integration is to be terminated. The second element of 
JLIMIT [ JLIMIT (2 ) ] denotes the allowable maximum number of radial 
stations when the mass flow rate exceeds the specified fraction of 
the estimated mass flow rate at X£ND(1). Each element of JLIMIT 
must be positive* odd, and no less than 5 but no greater than El. 
JLIHIT(l) should be less than or equal to JLIMIT(2). A default 
value of 13 is specified for JLIMIT(l), and a default value of 21 
is specified for JLIMIT{2). 


6. NAMELIST LIST4 

The parameters entered in namelist LIST4 specify the thermodynamic 
model, and the molecular transport properties. 

R A positive real variable denoting the gas constant, in either 

(ft-lbf)/(slug-R) or (J/kg-K). A default value of 1716.16116 
(ft-lbf)/(slug-R) is specified for R. 

GAMMA A positive real variable denoting the specific heat ratio. A de- 
fault value of 1.4 is specified for GAMMA. 

• If the influence of molecular transport is not to be included in the 
computation (KVISCY = 0), no other parameters must be specified in namelist 
LIST4. If the viscous, and thermal diffusion terms are to be included in the 
computation, the parameters presented in the following discussion must be 
specified. 

The dynamic viscosity is represented in the computer program by the 
Sutherland formula (2). 


It 

1 .5 

f T + B 
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T + B 
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In equation (1), y is the dynamic viscosity at the 
B is a constant. For air, B has the value 198.6 R 
is the viscosity at the reference temperature T„. 
must be specified in the program input by entering 
meters . 


absolute temperature T, and 
(1T0 K) . The parameter p 
The constants y p , T 0 , and°B 
the following, three para- 
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VIS0 A positive real variable denoting the reference viscosity y 0 in 
equation (1)-. The units of VIS0 are either (1 bf-sec/ft 2 } or 
(N-s/m 2 ). A default value of 3.5 x 10 -7 (lbf-sec/ft 2 ) is specified 
for VIS0 (this value' i's the dynamic viscosity of air at 492.0 R). 

T0 A positive real variable denoting the reference absolute temperature 

T 0 in equation (1). The units of T0 are either R or K. The specified 
value of T0 must correspond to the specified value of VIS0. A de- 
fault value of 492.0 -R is specified for T0. 

B A' positive' real variable denoting the constant B in equation (1). 

. The units for B are either R or K. A default value of 198.6 R is 
- ’ specified for B. 

Thermal conductivity is represented in the computer program by a quadratic 
curve fit written as a function of absolute -temperature. This quadratic poly- 
nomial is given by 


K = a l + a 2 ^ + a 3^ 2 


( 2 ) 


where k is the thermal conductivity, T is the absolute temperature, and the_ 
coefficients a. (i=l,2,3) are determined by fitting this expression to three 
data point sets. The curve fit coefficients are determined by entering the 
following two arrays. 

TDL A one-dimensi.onal real variable array consisting of three elements. 

Each element of TDL represents a temperature data point for the 
thermal conductivity of the working gas. The units of each element 
of TDL are either R or K. The successive elements of TDL must rep- • 
resent monotonically increasing, temperatures . The default values 
for the elements of TDL are: 

TDL ( 1 ) = 400.0 R 

TDL ( 2 ) = 1400.0 R 
TDL (3 } = 2400.0 R 

C0ND A one-dimensional real variable array consisting of three elements. 

Each element of C0ND denotes the working gas thermal conductivity at 
the absolute temperature specified by the corresponding element of 
TDL. The units of each element of C0ND are either (ft-slug)/(sec 3 -R) 
or (m-kg)/(s 3 -K). The default values for the elements of C0ND' are: 

C0ND ( 1 ) = 2.550963 x 10~ 3 (ft-slug)/(sec 3 -R) 

C0ND(2) = 7.566417 x 10~ 3 (ft-slug)/(sec 3 -R) 

C0ND(3) = 1.145772 x 10“ 2 (ft-slug)/(sec 3 -R) 
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The specified default values for the elements of C0ND correspond to 
the thermal conductivity values of air at the temperatures specified 
by the elements of TDL. 


7. NAMELIST LIST5 

The input parameters entered in namelist LIST5 specify the contours of 
the centerbody and the cowl. It is assumed that both the centerbody and the 
cowl are axisymmetric. The x-coordinate axis is the longitudinal axis of 
both the centerbody and the cowl (see Figure 1). The forebody tip must be 
located at x = 0.0. The axial station of the cowl lip must be specified at 
x > 0.0. 

For the purpose of geometry description, the axial (x) domain is divided 
into a number of intervals, as illustrated in Figure 5. The number of axial 
stations at which the centerbody geometry is specified is denoted by NCENT. 

The number of intervals on the centerbody is equal to (NCENT - 1). The number 
of axial stations at which the cowl geometry is specified is denoted by NC0WL. 
The number of intervals on the cowl is equal to (NC0WL - 1). 

In any interval, the centerbody or cowl radius may be specified either 
by tabular input, or b.y supplying the coefficients in a cubic polynomial writ- 
ten as a function of x . For the tabular input option., the body radius r(x) 
at axial position x in the ith interval (x| < x < Xi+l) is found by linear 
interpolation between point ,r-j } and point (xi+i ) . The local slope 

of the body for this interval in a given meridional plane is then given by 
the slope of the line segment joining these two points. Alternatively, 
employing the cubic polynomial 


r(x) = a. + b i (x - X j) + c^x - x^ 2 + d^x - x^.) 3 (x i < x < x- +1 ) (3) 


requires that the curve fit coefficients a-j, b-j , c-j, and d-f be supplied by the 
user. Since equation (3) is a cubic, slope and curvature can be matched at 
the junction point between two adjacent intervals employing this formulation. 
An option exists to employ the following cubic polynomial instead of equation 
(3). 


r(x) = a i + b^x + c^x 2 + d^.x 3 (x^ < x < x i+1 ) 


(4) 


When employing equation (3) or equation (4), the coefficients a-j, bj, c-j and- 
d-j must be specified for up to at most (NCENT - 1) and/or (NC0WL - 1) 
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intervals. The axial location must always be specified for NCENT and (if 
the internal flow field is to be computed) NC0WL axial stations, no matter 
which formulation is used. 

If the forebody is conical ahead of a certain axial station, the forebody/ 
centerbody geometry in this interval (1st interval) may be specified by 
entering the cone half-angle directly rather than by supplying the curve fit 
coefficients or entering the body radius by tabular input. 

The geometry description option for a- given interval is specified by the 
user and does not have to be the same for all intervals. For instance, the. 
forebody/centerbody contour may have a conical tip, then a quadratic or cubic 
variation with x, then a linear variation with x, then again a quadratic or 
cubic variation with x. A way to describe this contour would be to input 
the cone half-angle for the first interval, the cubic curve fit coefficients 
for the second interval, the body radii at the ends of the third interval, and 
the cubic curve fit coefficients for the fourth interval. A1 ternatively, the 
appropriate cubic curve fit coefficients could be supplied for each of the 
four intervals. It should be noted that in the input of the geometry data, 
radius, slope, and curvature should be made compatible between adjacent 
intervals. The boundary contours are specified by entering the following para- 
meters. 

KBASE An integer variable denoting whether equation (3) or equation (4) is 
to be employed when at least one interval of the centerbody or cowl 
geometry is - specified by either of these two cubic equations. If 
neither equation (3) nor equation (4) is used in the geometry de- 
scription, then KBASE does not have to be specified. If a cubic 
equation is to be used, then entering KBASE = 0 specifies that 
equation (3) will be employed. If KBASE = 1, then equation (4) will 
be used. A default value of 0 is specified for KBASE. The specified 
value of KBASE applies to both the forebody/centerbody contour and 
the cowl contour. 

NCENT A positive integer variable denoting the number of axial stations 
used in specifying the forebody/centerbody geometry. The number 
of intervals for the forebody/centerbody is equal to (NCENT - 1). 

The specified value for NCENT must be at least 2 and no greater 
than 100. A default value of 2 is specified for NCENT. 

KDCENT A one-dimensional integer variable array (dimensioned at 100) 

consisting of (NCENT - 1) elements. Each element of KDCENT specifies 
the forebody/centerbody geometry description option to be used in 
the corresponding interval. Specifying KDCENT(I) = 1 for the i th 
interval selects the option in which the forebody/centerbody radius 
is described by equation (3) or equation (4) (depending on the value 
of KBASE). Specifying KDCENT(I) = 2 for the ith interval selects 
the option in which the forebody/centerbody radius is specified by 
tabular input. If the forebody tip is conical, the geometry in 
the first forebody/centerbody interval may be described by entering 
the cone half-angle directly and specifying KDCENT(l) = 3. The 
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specified value for KDCENT(I) must be 1 or 2, except for the first 
interval (I = 1) in which case values of 1, 2, or 3 may be specified. 
Specifying other values than those allowed causes the program execu- 
tion to be aborted. The default value for KDCENT(l) is 3, while all 
other elements of KDCENT have no default values specified. 

XCENT A one-dimensional real variable array (dimensioned at 100} consist- 
ing of NCENT elements.. Each element of XCENT denotes the axial (x) 
position, in either feet or meters, of the beginning of a forebody/ 
centerbody interval [XCENT(NCENT) denotes the axial position of the 
end of the last interval]. The elements of XCENT must be nonnegative 
and montonically increasing. The default values for XCENT(l) and 
XCENT(2) are 1.0 ft and 3.5 ft, respectively. The remaining elements 
of XCENT do not have default values specified. 

RCENT A one-dimensional real variable array (dimensioned at 100) consist- 
ing of up to NCENT elements. Each element of RCENT specifies the 
forebody/ centerbody radius, in either feet or meters, at the axial 
location specified by the corresponding element of XCENT. If 
KDCENT(I) = 2, then RCENT(I) and RCENT (I + 1) must be specified. If 
KDCENT(I) =1 or 3, then RCENT(I) and RCENT(I + 1 ) do not have to be 
specified. Each element of RCENT must be nonnegative. No default 
values are specified for the elements of RCENT. 

ACENT One-dimensional real variable arrays (each array is dimensioned at 
BCENT 100), where each array consists of up to (NCENT - 1 ) elements. These 

CCENT arrays are used in conjunction with equation (3) or equation (4) for 

DCENT specification of the forebody/centerbody geometry. The elements of 
ACENT, BCENT, CCENT, .and DCENT specify the coefficients a-j , b-j , Cn, 
and d-j, respectively, in equation (3) or equation (4). If KDCENT(I) = 
1, then ACENT(I), BCENT(I) , CCENT(I), and DCENT ( I ) must be specified 
for the ith interval. If KDCENT(I) = 2 or 3, then ACENT(I), BCENT(I), 
CCENT(I), and DCENT(I) do not have to be specified for that interval. 
The units for the elements of ACENT are either feet or meters. The 
elements of BCENT are dimensionless. The units for the elements of 
CCENT are either (feet)" 1 or (meters) -1 . The units for the elements 
of DCENT are (feet) -2 or (meters)- 2 . No default values are specified 
for- the elements of ACENT, BCENT, CCENT, and DCENT. 

C0NE A real variable denoting the cone half-angle, ‘in degrees, of the 

forebody tip if it is conical. If KDCENT(l) = 3, then C0NE must be 
specified. If KDCENT(l) = 1 or 2, then C0NE does not have to be 
specified. A default value of 10.0 degrees is specified for C0NE. 

If only the external flow field about the forebody is to be computed 
[KCALL(l) = 1, KCALL(2) = 0, and KCALL(3) = 0], no further parameters must be 
specified in namelist LIST6. If the flow field in the annulus is to be 
determined, the following parameters must be entered. 
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NC0WL 


KDC0WL 


XC0WL 


RC0WL 


AC0WL 

BC0WL 

CC0WL 

DC0WL 


A positive integer variable denoting the number of axial stations 
used in specifying the cowl geometry. The number of intervals for 
the cowl is equal to (NC0WL - 1). The specified value for NC0WL 
must be at least 2 and no greater than 50. A default value of 2 
is specified for NC0WL. 

A one-dimensional integer variable array (dimensioned at 50) 
consisting of (NC0WL - 1) elements. Each element of KDC0WL speci- 
fies the cowl geometry description option to be used for the 
corresponding interval. Specifying KDC0WL ( I ) - 1 for the i th 
interval selects the option in which the cowT radius is described 
by equation (3) or equation (4) (depending on the value of KBASE). 
Specifying KDC0WL ( I ) = 2 for the ith interval selects the option in • 
which the cowl radius is specified by tabular input. The specified 
value for KDC0WL(I) must be either 1 or 2. Specifying other values 
than those allowed causes the program execution to be aborted. The 
specified default value for KDC0WL(1 ) is 2, while all other elements 
of KDC0WL have no default values specified. 

A one-dimensional real variable array (dimensioned at 50) consisting, 
of NC0WL elements. Each element of XC0WL specifies the axial (x) 
position, in either feet or meters, of the beginning, o.f a cowl 
interval [XC0WL(NC0WL) denotes the axial position of the end of the 
last interval]. Each element of XC0WL must be nonnegative and 
monotonical ly increasing. The default values for XC0WL(1), and 
XC0WL{2) are 2.0 ft and 3.5 ft, respectively. The remaining elements 
of XC0WL do not have default values specified. 

A one-dimensional real variable array (dimensioned at 50) consisting 
of up to NC0WL elements. Each element of RC0WL specifies the cow! 
radius, in either feet or meters, at the axial location specified by 
the corresponding element of XC0WL. If KD€0WL(I) = 2, then RC0WL(I) 
and RC0WL( I + 1) must be specified. If KDC0WL(I) = 1, then RC0WL(i) 
and RC0WL(I + 1) do not have to be specified. Each element of RC0WL 
must be nonnegative. The default values for RC0WL(1) and RC0WL(2) 
are both 0.70 ft. The remaining elements of RC0WL do not have 
default values specified. 

One-dimensional real variable arrays (each array is dimensioned at 
50), where each array consists of up to (NC0WL - 1) .elements. These 
arrays are used in conjunction with equation (3) or equation (4) 'for 
specification of the cowl geometry. The elements of AC0WL, BC0WL, 
CC0WL, and DC0WL specify the coefficients a-j, b-j, c-j, and d-j, 
respectively, in equation (3) or equation (4). If KDC0WL(I) = 1, 
then AC0WL(I), BC0WL(I), CC0WL(I), and DC0WL(I) must be specified 
for the ith interval. If KDC0WL(I) = 2, then AC0WL(I) J BC0WL(I), 
CC0WL(I), and DC0WL(I) do not have to be specified for that interval. 
The units of the elements of AC0WL are ei-ther feet or meters. The 
elements of BC0WL are dimensionless. The units for the elements of 
CC0WL are either (feet) -1 or (meters)- 1 . The units for the elements 
of DC0WL are either (feet)" 2 or (meters)- 2 . No default values are 
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specified for the elements of AC0WL, BC0WL, CC0WL, and DC0WL. 

DXTRAN A real variable denoting the centerbody translation from the design 
point position, or, equivalently, the amount the cowl has been trans- 
lated with respect to the centerbody. The units of DXTRAN are 
either feet or meters. Translation occurs solely in the x-direction. 
Moreover, the origin of the coordinate system is maintained at the 
forebody tip when translation occurs. A positive value for DXTRAN 
corresponds to a forward centerbody translation or a rearward cowl 
translation. A default value of 0.0 is specified for DXTRAN. 


8. NAMELIST LIST6 

The parameters entered in namelist LIST6 specify the various convergence 
tolerances and iteration limits used in the numerical integration. All para- 
meters in this namelist have specified default values. In general, the program 
is executed without changing the values of any of the parameters in this name- 
1 ist, 

SAFEIN A positive real variable denoting the ratio of the axial marching 
step taken to the axial marching step allowed by the Courant- 
' Friedrichs-Lewy (CFL) stability criterion. This variable is used 
to determine the axial position of both the first solution plane in 
the forebody flow field integration [KCALL(l) = 1] and the first 
solution plane in the internal flow field integration in which 
shock waves are not discretely fitted [KCALL(3) = 1]. Ensuing solu- 
tion planes for these integration options have their axial locations 
adjusted in accord with an internally computed value of SAFEIN. 

For the internal flow field integration option in which shock waves 
are discretely fitted [K.CALL(2) = 1], the axial position of each 
solution plane (except in the vicinity of a shock wave reflection) 
is controlled by the input value of SAFEIN. The specified value of 
SAFEIN must be positive, and must be less than 1.0 to satisfy the 
CFL stability criterion. The default and recommended value of 
SAFEIN is 0.975. 

CRIT A one-dimensional real variable array consisting of 16 elements. 

Each element of CRIT specifies a convergence tolerance or other 
parameter. The elements of CRIT have the following definitions 
and default values. 

CRIT(l) A positive real variable denoting the tolerance, in either feet or 
meters, used to determine if a user supplied initial-value plane 
data point is sufficiently close to a plane of symmetry when the 
data point is supposed to lie on the plane of symmetry. CRIT(l) is 
also used to determine if a user supplied initial -value plane data 
point is sufficiently close to the solid boundary when that point 
is supposed to lie on the solid boundary. A default value of ;0.1 ft 
(or 0.1 m) is specified for CRIT(l), 
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CRIT(2) 


CRIT(.3) 


CRIT(-4) 


CRIT(5) 


CRIT(6) 


CRIT(7) 


,CRIT(8) 


A positive real variable denoting the relative tolerance used in 
testing for a loss of significance in IBM library subroutine GELG 
(GELG is used to solve a system of simultaneous linear equations),. 

A default value of 10~ 7 is specified for CRIT(2). 

A positive real variable denoting the relative tolerance used in 
testing for convergence in the internal generation of the initial- 
value plane flow property field. A default value of 10“ 4 is 
specified for CRIT(3). 

A positive real variable denoting the relative tolerance used in 
testing for convergence of all three coordinates in the iterative 
scheme employed in computing a streamline-surface intersection or 
a bicharacteristic-surface intersection. A default value of 10~ 5 
is specified for CRIT(4). 

A positive real variable denoting the relative tolerance used in 
testing for the convergence of the five flow properties u, v, w» 

P, and p in subroutine S0LVE. A default value of 10“^ is specified 
for CRIT(5) . 

A positive real variable denoting the relative tolerance used in 
testing for the convergence of the static pressure in subroutine 
SH0CK. Convergence is attained in the local iteration loop if 


| P(2) - P*(2)|/P(2) < CRIT(6) 


where P(2) is the solution point pressure obtained from the local 
Hugoniot relations, and P*(2) is the pressure obtained from the 
wave surface compatibility relation. A default value of 10“ 4 is 
specified for CRIT(6). 

A positive real variable denoting the mass flow rate ratio at 
which the maximum number of radial stations allowed in the fore- 
body flow field computation is changed from JLIMIT(l) to JLIMIT(2). 
The mass flow rate ratio is the mass flow rate at a given fore- 
body flow field solution plane divided by the estimated mass flow 
rate at the axial station corresponding to XEND(l). A default 
value of 0.5 is specified for CRIT(7). 

A positive real variable used as a multiplier of the mass flow 
ratio which is employed in determining whether or not point 
addition is to be performed on a solution plane in the forebody 
flow field integration. The mass flow rate ratio is the mass flow 
rate at the solution plane just computed divided by that at the 
last solution plane where point addition or deletion was performed. 
A default value of 1.0 is specified for CRIT(8). 
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CRIT(9) 


CRIT(IO) 


CRIT (11) 


CRIT (12) 


A positive real variable denoting the relative- tolerance used in 
routine LINK21 for determining when the angle a calculated in the 
global correction for the bow shock wave points has converged. A 
default value of 10 -4 is specified for CRIT('9). 

A positive real variable used in routine LINK21 for determining- if 
a sufficient number of shock wave solution points have converged in 
global correction. Convergence is attained when 


M/IST0P > CRIT ( 10) 


where M is the number of shock wave solution points which have 
converged in global correction, and IST0P is the number of cir- 
cumferential stations in the computed sector. A default value of 
0.8 is specified for CRIT (10)- 

A positive real variable used in subroutine BSH0CK for determining 
if the velocity component downstream of the shock wave and normal 
to the surface of the solid boundary has converged to within a 
specified tolerance of zero. Convergence is attained when 


fV N | < CRIT(ll) 


where is the velocity component normal to the solid boundary. 

A default value of 10 -6 ft/sec (or 10 -6 m/s) is. specified for 
CRIT(l.l). 

A positive real variable used in routine LINK22 for determining if 
another solution plane is to be inserted between the last solution 
plane and the intersection of the incident internal shock wave with 
the solid boundary. Another solution plane is inserted if 


Ax/Ax cfl > CRIT{12) 


where Ax is the axial (x) distance between the last computed plane 
and the nearest point on the space curve defined by the intersec- 
tion of the incident internal shock wave with solid boundary, and 
AXcpi is the axial step allowed by the Courant-Friedrichs-Lewy 
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CRIT-(1 3) 
CRIT( 14) 


CRITO 5) 


CRIT(16) 


(CFL) stability criterion, A default value of 0.2 is specified for 
CRITO 2) . 

Not presently employed. 

A positive real variable used in subroutine PENTRE for determining 
if a streamline-shock wave intersection point is sufficiently close 
to the current solution plane, so that an interior point unit process 
on the downstream side of the shock wave is not performed. Instead, 
a streamline projection onto the solution plane and subsequent flow 
property interpolation in this plane is performed. The application 
of the interior point unit process is not performed if 


(x s - x int^ /Ax CFL < CRIT -( 14 ) 


where x is the axial position of the solution plane, -x. . is the 
axial location of the streamline-shock wave intersection-point, and 
Axcfl is the axial marching step allowed by the Courant-Friedrichs- 
Lewy stability criterion. A default value of 0.4 is specified for 
CRIT( 1 4} . 

A positive real variable used in subroutine STRSHK for determining 
if convergence has been obtained in calculating the intersection 
point .of a body streamline with the space curve defined by the 
intersection of the incident internal shock wave with a solid bound- 
ary. Convergence is attained when 


l 0 i+ l “ e il < CRITO 5) 


where 6. is the polar angle of the intersection point on the Tth 
iteration, and 0j + -j is the polar angle on the (i+l)tf^ iteration. 

A default value of 10 -4 radians is specified for CRIT (15). 

A positive real variable used in subroutine INTSCT for determining if 
convergence has been obtained in calculating the intersection point 
of a bicharacteristic with either a solid boundary or a shock wavfes 
or the intersection point of a streamline with a shock wave. Conver- 
gence is attained when 


- R s | < CRIT(1 6) 
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where R-j is the radius of the intersection point obtained by inte- 
grating the equation for a streamline or bicharacteristic, and R s 
is the intersection point radius obtained from the shock wave or 
boundary surface formulations. A default value of 10" G ft (or 
10" 6 m) is specified for CRIT(16). 

ITEND A one-dimensional- integer variable array consisting of 6 elements. 

Each element of ITEND specifies a. limit to the number of iterations 
permissible in a given iteration loop. The elements of ITEND have 
the following definitions and default values. 

ITEND(l) A positive integer variable denoting the maximum number of inner 
iterations permissible in determining, the intersection coordinates 
of either a streamline with a surface, or a bicharacteristic with a 
surface. ITEND(l) is used in conjunction- with CRIT(4). A default 
value of 10 is specified for ITEND(l), 

ITEND(2) -A positive integer variable denoting the maximum number of outer 
iterations permissible in obtaining the five flow properties u, v, 
w, P, and p in all unit processes except the shock wave-solid boundary 
point unit process. A default value of 10 is specified for ITEND(2). 

ITEND(3} Not presently employed. 

ITEND(4) A positive integer variable denoting the maximum number of iterations 
permissible in the relaxation of the velocity component normal to 
the solid boundary and downstream of the reflected (cowl lip) shock 
wave in the shock wave-solid boundary point unit process (subroutine 
BSH0CK). A default value of 20 is specified for ITEND(4). 

ITEND(5) Not presently employed. 

ITEND(6) A positive integer variable denoting the maximum number of permissible 
subiterations in determining the intersection point of a line segment 
with a given three-dimensional surface (subroutine INTSCT). A default 
value of 10 is specified for ITEND(6). 


9. NAMELIST LIST7 

The parameters entered in namelist LIST7 specify if debug output is to be 
printed. 

KDUMP A one-dimensional integer variable array consisting of four elements. 

Each element of KDUMP specifies whether or not a particular unit proc- 
ess is to have debug output printed. Specifying KDUMP(I) = 1 (1=1 to 
4) activates the debug output option for the corresponding unit 
process. Specifying KDUMP(I) = 0 causes no debug output to be printed 
for the corresponding unit process. The elements of KDUMP activate 
the debug output option for the following unit processes and have the 
following default values. 
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KDUMP(I) 

Activates Debug Output for 

Default Value 

KDUMP(l) 

interior point scheme 

0 

KDUMP(2) 

solid body point scheme 

0 

KDUMP( 3 ) 

field-shock wave point 
scheme 

0 

KDUMP{4) 

solid body-shock wave 
point scheme 

0 


ISTART A positive integer variable denoting the solution plane number at 
which debug output is to be initiated. A default value of 1 is 
specified for ISTART. 

10. FORMATTED READ OF THE INITIAL-VALUE PLANE FLOW PROPERTY FIELD 

The user supplied initial -value plane flow property field is entered by 
a formatted read of file ITAPE after all seven namelists have been input. To. 
enter the initial -val ue plane by tabular input, KIVS = 0 must be specified in 
namelist LIST2. The default value for ITAPE is 5(the input file). The 
initial -value plane may be read from TAPE 3 by specifying ITAPE = 3 in name- 
list LIST2. TAPE 3 is linked to the dummy file IVS in the PR0GRAM card. 

The index limits for the initial-value plane, IST0P and JMAXI, are 
entered in namelist LIST3. The initial-value plane point networks for the 
four flow symmetry options are illustrated in Figure 3. 

The initial -value data are entered by the formatted read statement 

READ (ITAPE, 260) ((Y(I ,J) ,Z(I ,J) ,U( I,J) ,V(I,J) ,W(I,J) ,P(I,J) ,RJ9(I,J) , 

J=1 , JMAXI), 1=1, IST0P) 

with the format (4E20.13/3E20.13) . The parameters in the formatted read state 
ment have the following definitions (see Figure 3). 

I An integer denoting the circumferential index of the data point. 

J An integer denoting the radial index of the data point. 

Y A two-dimensional real variable array, each element of which 

denotes the y-position, in either feet or meters, of a data 
point. The Y array is dimensioned at (40,21). No default 
values are specified for the elements of Y. 

Z A two-dimensional real variable array, each element of which 

denotes the z-position, in either feet or meters, of a data 
point. The Z array is dimensioned at (40,21). No default 
values are specified for the elements of Z. 

U A two-dimensional real variable array, each element of which 

denotes the x- component of velocity, in either (ft/sec) or 
(m/s), at a data point. The U array is dimensioned at (40,21), 

No default values are specified for the elements of U. 
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V A two-dimensional real variable array, each element of which denotes 

the y-component of velocity, in either (ft/sec) or (m/s), at a 
data point. The V array is dimensioned at (40,21). No default 
values are specified for the elements of V. 

W A two-dimensional real variable array, each element of which denotes 

the z-component of the velocity, in either (ft/sec) or (m/s), at 
a data point. The W array is dimensioned at (40,21). No default 
values are specified for the elements of W. 

P A two-dimensional real variable array, each element of which denotes 

the pressure, in either (lbf/ft 2 ) or (N/m 2 ), at a data point. The 
P array is dimensioned at (40,21). No default values are specified 
for the elements of P. 

R0 A two-dimensional real variable array, each element of which denotes 

the density, in either (slug/ft 3 ) or (kg/m 3 ), at a data point. The 
R0 array is dimensioned at (40,21). No default values are specified ■ 
for the elements of R0. 

In all cases, the initial-value plane data points with J = 1 must lie 
(to a close approximation) on the surface of the forebody/ centerbody. The 
initial-value plane data points with J = JMAXI correspond to the downstream 
bow shock wave points. If only the internal flow is to be computed, the 
initial-value plane is located at the cowl lip axial station. It is sufficient 
to specify the flow property field, in this case, to a point just outside of 
the cowl lip. If the bow shock wave radius is less than that of the cowl lip, 
the execution is aborted. 

For the case of no planes of flow symmetry (KSYM = 0), the computed 
sector is the entire solution plane corresponding to IST0P circumferential 
stations and JMAXI radial stations [see Figure 3(a)]. For the case of one 
plane of flow symmetry (KSYM =1), the computed sector is the half-plane 
bounded by the y-axis and containing the +z-axis [see Figure 3(b)], In this 
case, the data points with I = 1 must lie on the +y-axis, and the data points 
with I = IST0P must lie on the -y-axis. For the case of two planes of 'flow 
symmetry (KSYM =2),. the computed sector is the quadrant bounded by the +y-axis 
and the +z-axis [see Figure 3(c)]. In this case, the data points with I = 1 
must lie on the +y-axis, and the data points with I = IST0P must lie on the 
+z-axis. For the axi symmetric flow case (KSYM = 3), the computed sector is 
limited to the single circumferential station lying on the +y-axis [see 
Figure 3(d)]. In this case, the data points with I = 1 must lie on the +y-axis. 

If the forebody flow field is not being calculated [KCALL(l) = 0], or if 
it is being computed and the bow shock wave is conical (KC0N = 1), no other 
input parameters have to be entered. If, however, the forebody flow field is 
to be calculated [KCALL(l) = 1] and the bow shock wave is not conical (KC0N = 0), 
then the angle subtended by the bow shock wave and the x-axis in the meridional 
plane defined by the shock wave point for each shock wave point in the computed 
sector must be entered by the formatted read statement 
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READ (ITAPE,270) (BETA(I) ,1=1, IST0P) 


with the format (E20.13). The parameters in .the formatted read statement have 
the following definitions (see Figure 3). 

I An integer denoting the circumferential index of the initial-value 

plane' shock wave point. 

BETA A one-dimensional real variable array (dimensioned at 4Q) ? each 
element of which denotes the angle, in radians, subtended.by the 
bow shock wave and the x-axis in the meridional plane defined by 
the corresponding initial -value plane downstream shock wave point. 
Each element of BETA must be positive. No default values are 
specified for the elements of BETA. 
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SECTION V 


OUTPUT INTERPRETATION 


The initial portion of the computer output comprises preliminary infor- 
mation. This preliminary output consists of information identifying the 
problem being considered, the specified computation options, the flow symmetry 
option, the thermodynamic model and the molecular transport properties, the 
vehicle orientation and the free-stream conditions, certain index parameters, 
the contours of the centerbody and the cowl, and the convergence tolerances 
and the iteration limits. The initial -value plane is then printed. Alter- 
natively, if a program restart is specified, the last solution piane written 
on the restart file is printed. Each solution plane is then printed in a 
format similar to the initial-value plane printout. Additionally, the redis- 
tributed data plane at the cowl lip axial station is printed if the internal 
flow integration option is specified. Moreover, for the internal flow field 
computation, the solution points are printed which lie along the space curves 
defined by the intersection of the internal shock wave with the solid boun- 
daries. 

The output parameters listed on the computer printout are defined below. 

I circumferential index of the solution point 

J radial index of the solution point 

X axial position of the solution plane or the solution point 

Y y-position 

Z z-position 

M Mach number 

Q velocity magnitude 

P pressure 

R0 density 

T absolute temperature 

U x-component of velocity 

V y-component of velocity 

W z-component of velocity 

PT stagnation pressure 

TT stagnation temperature 

ITG number of global corrector applications 

ITL number of local iterations 
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Streamline solution points have both I and J indices which are numbers. 
An upstream shock wave solution point is denoted by a numerical I index and 
the J index is U. A downstream shock wave solution point is denoted by a 
numerical I index and the J index is D. 

For the external flow field about the forebody , the body streamline 
solution points are denoted by J = 1 . The outer bound to the computational 
flow regime is defined by the locus of downstream shock wave solution points, 
points with J = D. Since periodic point addition and deletion are performed 
in the external flow field integration, continuous streamlines throughout 
the computational flow regime are not available. Inserted solution, points 
are noted by ITG = ITL = 0. 

For the continuous internal flow field integration option, the body 
streamline points on the surface of the centerbody are denoted by J = T, and 
the body streamline points on the surface of the cowl are denoted J, = J INLET. 
The solution is found on the continuous streamlines which pass through 
the redistributed points on the solution plane at the cowl lip axial' station. 

For the internal flow field integration option in which shock waves are 
discretely fitted, the body streamline points on the surface of the- center- 
body are denoted by J = 1. The body streamline points on the surface of the. 
cowl are denoted by J = (JINLET - 2). The shock wave solution points float 
in the storage arrays as the internal shock wave travels between the center- 
body and the cowl on successive solution planes. The streamline points 
between the upstream side solid boundary (either centerbody or cowl) and' the: 
upstream shock wave points on a given solution plane' lie in the upstream’ 
flow field sector on that solution plane. In a like manner, streamline 
points which lie between the other solid boundary and the downstream - shock, 
wave points on a given solution plane Tie in the downstream flow field sector 
on that solution plane. A reversal of the upstream and downstream sectors 
occurs at an internal shock wave-solid boundary intersection. It should: be 
noted that continuous streamlines are followed in the internal flow field 
integration. 

The intersection of the incident internal shock wave with a solid boun- 
dary at a shock wave-solid boundary intersection defines a space curve. The 
solution is found on both the upstream and downstream sides- of both the in- 
cident and reflected shock waves at points on this space curve-. 

At the end of a solution plane printout, the Courant number and the- 
x-step regulation parameters are printed. The Courant number is the ratio of 
the axial step taken to the axial step allowed by the Courant-Friedrichs- 
Lewy stability criterion (based on immediate neighbors in the interpolation- 
fit point stencils). The Courant number listed is that used in obtaining 
the solution plane that was just printed. Likewise, the x-step regulation 
parameters refer to the solution plane that was just printed. 
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SECTION VI 


SAMPLE CASES . 


1 . INTRODUCTION 

Five sample cases are presented in this section to illustrate the applica- 
tion of the computer program for calculating the flow field in supersonic mixed- 
compression aircraft inlets. For each of the five sample cases, a discussion 
of the problem is given, the required input data are presented, and selected • 
portions of the computer output are listed. The input parameter discussions 
follow the order in which the input parameters are presented in Section IV. 

Sample Case No. 1 considers the computation of both the external flow field 
and the internal flow field for a simplified geometry mixed-compression inlet 
at zero angle of attack. Sample Case No. 2 considers the computation of the 
external flow field about the forebody of the inlet of Sample Case No. 1 when 
the angle of attack is 2.5 degrees. Sample Case No. 3 considers the calculation 
of the ax i symmetric internal flow field in a simplified geometry annulus using 
the program options in which internal shock waves are not discretely fitted 
but the molecular transport terms are included. Sample Case No. 4 considers 
the computation of the axisymmetric internal flow field in the Mach 3.5 inlet, 
documented in Reference (4), at' zero angle of attack. Sample Case No. 5 con- 
siders the computation of the internal flow field for the Mach 3.5 inlet for the 
off-design Mach number of 2.5, nonzero centerbody translation, and an angle of 
attack of 3.0 degrees. 


2. SAMPLE CASE NO. 1 

This sample case is concerned with the computation of both the external 
flow field and the internal flow field in a simplified geometry supersonic 
mixed-compression aircraft inlet. The inlet geometry is axisymmetric and the 
specified angle of attack is zero, hence the flow field is axisymmetric. This 
sample case represents the problem being considered when all program input 
parameters retain their default values. 

The first card of the data deck for Sample Case No. 1, which is presented 
in Figure 6, is the title card. All input parameters in namelist LIST! re- 
tain their default values. Thus, KUNIT = 1 and English units are used in the 
computation. The default values for KCALL(l), KCALL(2) , and KCALL{3) are 1, 

1, and 0, respecti vely. Consequently, the specified computation options are 
the forebody flow field integration option and the internal flow field inte- 
gration option in which shock waves are discretely fitted. The default values 
for XEND(l) and XEND(2) are 2.0 ft and 3.5 ft, respectively. Hence, the fore- 
body flow field integration terminates at x = 2.0 ft, and the internal flow 
field integration terminates at x = 3.5 ft. Note that since the internal flow 
field integration option in which shock waves are not discretely fitted is not 
used [ KCALL ( 3 ) = 0], XEND(3) does not have to be specified. The default value 
of RCAVG = 0.8 ft is used for estimating the mass flow rate downstream of the 
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Figure 6. Data deck for’ Sample Case No. 1. 



bow shock wave at the cowl lip axial station. The default value of KVISCY is 
0, consequently the molecular tranport terms are not included in the computa- 
tion. KSYM retains its default value of 3, hence the axisymmetric flow option 
is employed. The default value of KSGL0B = 1 specifies that global correction 
is to be performed in obtaining the solution for the bow shock wave points. 

KPRINT retains its default value of 1, consequently all solution points are 
printed. Since IPRSTP and KSTART are both 0, the execution is not terminated 
at a specified solution plane, nor are any restart file operations performed. 

All input parameters in namelist LIST2 retain their default values. Con- 
sequently, the freerstream Mach number MFS, the free-stream pressure PFS, and 
the free-stream density R0FS have values of 3.0, 242.2 (lbf/ft 2 ), and 0.0003622 
(slug/ft 3 ), respectively. PITCH and YAW retain their default values of 0.0, 
so the pitch and yaw angles are both 0.0. The axial position of the initial- 
value plane, specified by the default value of XI, is 1.0 ft. Since KIVS and 
KC0N are both 1, the initial -value plane is internally generated and the bow 
shock wave is assumed to be conical (the forebody is conical). The parameter 
ITAPE is not employed since the initial -value plane is generated internally 
(KIVS =1). 

All input parameters in namelist LIST3 retain their default values. 
Consequently, IST0P, JMAXI , and JINLET have values of 1, 11, and 11, respec- 
tively, so that 1 circumferential station is employed, 11 radial stations 
on the initial-value plane are specified, and 11 radial stations in the- inter- 
nal flow field computation are specified. JLIMIT(l) and 0LIMIT(2) retain 
their default values of 13 and 21 radial stations, respectively. 

All input parameters in namelist LIST4 retain their default values. 

Thus, the specific heat ratio and gas constant, specified by GAMMA and R, 
respectively, have values of 1.4 and 1716.16116 (ft-lbf )/(slug-R) , respectively. 
Since the molecular transport terms' are not included in the computation 
(KVISCY = 0), the input parameters VIS0, T0, B, TDL, and C0ND are not employed. 

All input parameters in namelist LIST5 retain their default values. The 
default inlet geometry has a conical forebody/ centerbody with a cone half-angle 
of 10.0 degrees. The forebody tip is located at x = 0.0 ft, and the centerbody 
geometry is specified to x = 3.5 ft. Thus, NCENT = 2, KDCENT(l) =3, 

XCENT(l) = 1.0, XCENT (2) = 3.5, and C0NE * 10.0. The cowl lip axial station is 
located at x = 2.0 ft and the cowl geometry is specified to x = 3.5 ft. The 
cowl radius, entered by tabular input, is constant at the value of 0.7 ft. 

Thus, NC0WL = 2, KDC0WL(1) = 2, XC0WL(1) = 2.0, XC0WL(2) = 3.5, RC0WL(l) = 0.7, 
and RC0WL(2) = 0.7. The remaining parameters in namelist LIST5 (RCENT, ACENT 
to DCENT, AC0WL to DC0WL, KBASE, and DXTRAN) are not employed. 

All convergence tolerances and iteration limits retain their default values 
in namelist LIST6. 

No debug output is to be printed, hence all input parameters in namelists 
LIST7 retain their default values. 

Selected portions of the computer output for this sample case are presented 
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in Figure 7. The first three pages of the program output present the job title, 
the specified computation options, the flow symmetry option, the thermodynamic 
model, the vehicle orientation and free-stream 'data, the type of initial -value 
plane, certain index parameters, the centerbody and the cowl contours, and the 
various convergence tolerances and iteration limits. The fourth page presents 
the internally generated initial -value plane flow property field. 'Pages five 
to seven present the first three solution planes on the forebody (planes No. 1 
to No. 3), and page 8 presents the last forebody solution plane (plane No. 16), 
which is located at the cowl entrance. Page nine presents the redistributed 
flow property field at the cowl entrance, and page ten presents the first 
internal flow field solution plane (plane No. 17). Pages 11 to 13 present the 
internal flow field solution plane just upstream of the first shock wave inters 
section with the centerbody (plane No. 26), the interplanar results, and the 
solution plane just downstream of the shock wave-centerbody intersection 
(plane No. 27). Pages V4 to 16 present the corresponding results for the 
first shock wave intersection with the cowl. The last page presents the last 
solution plane (plane No. 78) before the flow becomes subsonic. 

Sample Case No. 1 required approximately 95 seconds of central processor 
time on the CDC-6500 computer (using the Purdue University modified 6OOO-SC0PE 
compiler). 
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THE ANALYSIS OF STEADY THREE-DIMENSIONAL FLOW IN SUPERSONIC MIXED-COMPRESSION AIRCRAFT INLETS 


ABSTRACT 

THE FLOW FIELD IN A SUPERSONIC MIXED-COMPRESSION AIRCRAFT INLET IS COMPUTED USING THE METHOD OF CHARACTERISTICS 
FOR STEADY THREE-DIMENSIONAL FLOW. THE BOw SHOCK WAVE AND REFLECTED INTERNAL SHOCK WAVE SYSTEMS ARE' COMPUTED 

USING A DISCRETE SHOCK-FI TTING PROCEDURE. THE PROGRAM HAS THE CAPABILITY TO INCLUDE THL INFLUENCE OF MOLECULAR 

TRANSPORT ON THE SOLUTION BY TREATING THESE EFFECTS AS CORRECTION TERMS IN THE CHaRACTERISIICS SCHEME. 

THIS PROGRAM WAS DEVELOPED AT THE PURDUE UNIVERSITY THERMAL SCIENCES AND PROPULSION CENTER BY J. -VADYAK UNDER 
N.A.S.A, GRANT No, NGR- 15-005*191 FOR THE N.A.S.a. LEWIS RESEARCH CENY ER i CLEVELAND, OHIO, THE PRINCIPAL 
INVESTIGATOR WAS J.D, HOFFMAN ANO THE N.A.S.A. TECHNICAL DIRECTOR WAS A. BISHOP. 


o 

'S'Q 


UOB TITLE 

SAMPLE CASE NO. 1 . 


SPECIFIED COMPUTATION OPTIONS 
l.J FOREBODY FLOW FIELD 

8.) INTERNAL FLOW FIELO WITH SHOCK WAVE SYSTEM 

FLOW SYMMETRY 

AXISYMMETRIC FLOW 

GLOBAL CORRECTION 

GLOBAL CORRECTION IS PERFORMED ON THE BOW SHOCK WAVE POINTS 
THERMODYNAMIC MODEL 

A THERMALLY AND CALORICALLY PERFECT GAS IS SPECIFIED WITH 

SPECIFIC HEAT RATio=l.**o6od GAS C0NSTANT= 1 .716161E+03 ( FT-LBF /SLUG-DE6 Rj 

VISCOSITY AND THERMAL CONDUCTIVITY TRANSPORT TERM* 

VISCOUS AND THERMAL DIFFUSION TERMS ARE NOT INCLUDED IN THE COMPUTATION - INVISCIO ANO ADIABATIC FLOW IS ASSUMED 


ORIENTATION AND FREE STREAM DATA 

ORIENTATION • PITCHs 0 . 00000 ( DEGREES ) YAWs 0,000 00 {DEGREES} 

Figure 7. Selected output from Sample Case No. 1. 
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FREE STREAM DATA • MACH NO,= 3< 00000 PRESSURES £ ,R£200 DE+02 ( LBF7FT**2 ) DENSITY:: 

TEMPERATURES 3.896437E«-Q2<DEG Rl SONIC SPEEDS 9 .675577E+02 ( FT/SEC ) 

X-VEEOCITY— 2.9O2673E+0 3 (FT/sEC) Y-VELOCITYs 0» ' (FT/SEC) 2»VELOClTTc 0. 


INITIAL VALUE SURFACE 

an INTERNALLY GENERATED INITIAL VALUE SURFACE IS SPECIFIED AS BEING LOCATED AT Xs 1 .0 00 OOOt+O 0 ( FT) 
A CONICAL BOH SHOCK HAVE IS SPECIFIED - ThE INTERNALLY 6ENERATEO SHOCK WAVE ANGLES ARE 
BETA! 1>3 3.709GS7E-O1IRADIANS) 

INOEX PARAMETERS 

ISTOPs 1 JMAX=A0 JMAXJsll 

ULZMIT 11 ) s 13 ULIMIT ( 2>s21 

JINLET=U 


INTEGRATION TERMINATION POINTS 

FOREBOdY FLOW FIELD INTEGRATION TERMINATES AT X= 2,OOOOODE+00<FT) 
INTERNAL FLOW FIELD INTEGRATION TERMINATES AT X* 3.50000CE+00 I FT > 
FOR FOREBODY FLOW • RCAVG= 8. OOOOOOE-Ol < FY ) 
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Figure 7. Continued 



CONVERGENCE TOLERANCES i ITERATION LIMITS. AND OTHER PARAMETERS 

CONVERGENCE TOLERANCES AND OTHER PARAMETERS 

CRIT ( Its l.OQOOOOE-Ol 
CRIT( 21 = l.OOBOOOE-07 
CRIT ( 31= l.OUCOOOE-OR 
CRIT) 41 = I.0B00BBE-O5 
CRIT < 51 = 1.000000E-04 
CRIT) 6 1 c 1.000000E-04 
CRIT( 71= 5,OBOOOOE-Ol 
CRIT) 61= l,OOOUOOE.+aO 
CRIT ( 91= i,O0OO0OE-04 
CRIT) 10 ) = 6.000000E-01 
CRIT ( 11 1= 1.000000E-06 
CRIT ( 12 1= 2, OOOOOOE-O! 

CRIT ( II 1= 1.000000E-05 
CRIT 1 14 1 = 4,OUOOOOE-Ol 
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3. SAMPLE CASE NO. 2 


Sample Case No. 2 is concerned with the computation of the external flow 
field about the forebody of the inlet considered in Sample Case No. 1 when the 
angle of attack is 2.5 degrees. The initial-value plane flow property field 
is generated internally in the program by using the approximate technique 
described in Section II. 

The first card of the data deck for this sample case, illustrated in 
Figure 8, is the title card. English units are' used, so KUNIT retains its de- 
fault value of 1 in namelist LIST!. Since only the forebody flow field is to 
be computed, KCALL (2 ) is specified as 0 in namelist LIST1 , while KCALL(l) and 
KCALL(3) retain their default values of 1 and 0, respectively. The forebody 
flow field integration termination point, denoted by XEND(l), is 2.0 ft, the 
default value. Since one plane of flow symmetry exists, KSYM is specified as 
1. The molecular transport terms are not to be included in the computation, 
and global correction is to be performed on the bow shock wave solution points 
Hence, KVISCY and KSGL0B retain their default values of 0 and 1, respectively. 
All remaining parameters in namelist LIST! (RCAVG, KPRINT, IPRSTP, and KSTART) 
are left at their default values. 

In namelist LIST2, MFS, PFS, and R0FS are left at their default values of 
3.0, 242.2 (lbf/ft 2 ), and 0.0003622 (slug/ft 3 ), respectively. Since the .angle 
of attack is 2.5 degrees, PITCH is specified as 2.5. The remaining parameters 
in namelist LIST2 (YAW, XI, KIVS, KC0N; and ITAPE) retain their default values 
Consequently, the forebody flow field integration starts at x = 1.0 ft, the 
initial -value plane is internally generated, and a conical bow shock wave is 
assumed . 

In namelist LIST3, all input parameters retain their default values 
except for IST0P. The number of circumferential stations is selected to be 
21 , hence IST0P = 21 . 

All input parameters in namelist LIST4, which specifies the thermodynamic 
model, retain their default values. 

In namelist LISTS, only the forebody geometry must be specified [up to 
XEND(l)] since the internal flow field is not being computed. The default 
forebody/centerbody geometry is again used, hence NCENT ~ 2, KDCENT(l) - 3, 
XCENT(l) =1.0, XCENT(2) = 3.5, and C0NE = 10.0. No other parameters are used 
in namelist LIST5 for this sample case. 

All convergence tolerances and iteration limits specified in namelist 
LIST6 retain their default values. 

No debug output is to be printed, so all input parameters in namelist 
LIST7 are left at their default values. 

Selected portions of the computer output for Sample Case No. 2 are pre- 
sented in Figure 9. The preliminary information, the initial -value plane, 
and the first and last solution planes are presented. The execution of Sample 
Case No. 2 required 411 seconds of central processor time on the CDC-6500 
computer. 
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Figure 8. Data deck for Sample Case No. 2. 
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Figure 9. Continued 
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4. SAMPLE CASE NO. 3 


This sample case is concerned with the computation of the axisymmetric 
internal flow field for a simplified geometry inlet using the program option 
in which internal shock waves are not discretely fitted. For this sample case, 
the molecular transport terms are included in the computation of the internal 
flow field. 

The data deck for this sample case is presented in Figure 10. The first 
card' of the data deck is the title card. English units are used in the com- 
putation, hence KUNIT retains its default value of 1 in namelist LIST!. Since 
only the internal flow field integration option in which shock waves are not 
discretely fitted is to be used, KCALL(l) = 0, KCALL(2) = o, and KCALL(3) = 1 
are specified in namelist LIST1 . ■ XEND(3), which denotes the termination point 
of the internal flow field integration, is left equal to its default value of 
3.5 ft. To illustrate the use of the flow symmetry option in which two planes 
of symmetry exist (even though the flow field is axisymmetric), KSYM = 2 is 
specified in namelist LIST!. Since the transport terms are to be included in 
the computation, KVISCY = 1 is specified in namelist LIST!. The remaining in- 
put parameters in namelist LIST! are retained at their default values. Note 
that RCAVG and KSGL0B are not employed in this sample case. 

All parameters in namelist LIST2 are- left at their default values, except 
for XI, the axial location of the initial -value plane, which is specified as 
2.0 ft (the axial position of t|ie cowl lip). The initial -value plane flow 
property field is internally generated at x = 2.0 ft (a conical forebody is 
assumed) . 

In namelist LIST3, all input parameters retain their default values ex- 
cept for IST0P. The desired number of circumferential stations in the com- 
puted sector (quadrant) is 11, hence IST0P = 11. Both JMAXI and JINLET are 
11, hence the number of radial stations on both the initial -value plane and on 
each of the internal flow field solution planes is 11. Note that JLIMIT is 
not employed in the computation of an internal flow field. 

All parameters in namelist LIST4 are left at their default values. The 
default values of GAMMA = 1 .4 and R = 1716.16116 (ft-1 bf)/(slug-R) specify the 
thermodynamic model. The default values of VIS0 = 3.5 x 10" 7 (1 bf-sec/ft 2 ) , 

T0 = 492.0 R, and B = 198.6 R specify the parameters in Sutherland's law. The 
tabular values of TDL and C0ND presented in Section IV. 6 specify the thermal 
conductivity. 

The centerbody geometry for this inlet is identical to that specified in 
Sample Cases No. 1 and No. 2. Hence-, the default values of NCENT = 2, 

KDCENT(l) = 3, XCENT(l) = 1.0, XCENT(2) = 3.5, and C0NE = 10.0 are retained. 

The cowl geometry in this inlet is selected so that little or no flow turning 
occurs at the cowl lip. This requirement is met by specifying NC0WL = 2, 
KDC0WLO) = 2, XC0WL ( 1 ) = 2.0, XC0WL(2) = 3.5 , RC0WL(1 )=0. 7959, and RC0WL(-2) = 
0.8874 in namelist LIST5. Note that except for RC0WL(1) and RC0WL{2), these 
values are the default values. No other parameters are used in this namelist. 
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Figure 10. Data deck for Sample Case No. 3. 
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All convergence tolerances and iteration limits specified in namelist 
LIST6 retain their default values. 

No debug output is to be printed, hence all input parameters in namelist 
LIST7 are left equal to their default values. 

Selected portions of the computer output for Sample Case No. 3 are- pre- 
sented in Figure 11. The central processor time required to execute Sample 
Case No. 3 was 433 seconds on the CDC-6500 computer. 



The analysis OF STEADY THREE-DIMENSIONAL flow in SUPERSONIC MIXED “COMPRESS'! on aircraft inlets 
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abstract 

the flow fielo IN A SUPERSONIC MIXED-COMPRESSION aircraft inlet is computed using the method of CHARACTERISTICS 
for steady three-dimensional flow. the bow shock wave and reflected INTERNAL shock wave systems are computed 
using a discrete shock-fitting procedure. THE program has the capability to include the influence of molecular 

TRANSPORT ON THE SOLUTION BY TREATING THESE EFFECTS as CORRECTION terms in the CHARACTERISTICS scheme. 

THIS PROGRAM WAS DEVELOPED AT THE PURDUE UNIVERSITY THERMAL SCIENCES AND PROPULSION CENTER BY J. VaDYaK UNDER 
N.A.S.A. GRANT No. NGR-15-005-19I FOR THE N.A.S.A* LEWIS RESEARCH CENTER. CLEVELaNO. OHIO. THE PRINCIPAL 

INVESTIGATOR WAS J.O. HOFFMAN AND THE N.A.S.A. TECHNICAL DIRECTOR WAS A. BISHOP. 


JOB TITLE 

Sample case no. 3 


SPECIFIED COMPUTATION OPTIONS 

i.» internal flow fielo without shock wave system 


FLOW SYMMETRY 

TWO PLANES OF SYMMETRY - COMPUTED SECTOR IS THE OUAPRANT 80UN0ED BY THE *Y-AXIS AND THE +2-AXIS 


THERMODYNAMIC model 

A THERMALLY AND CALORICALLY PERFECT GAS IS SPECIFIED WITH 

SPECIFIC HEAT RATIO=I,«tOOOO GAS CONSTANTS 1 .71616IE+03 ( FT-LBF/SLUG-OEG R) 

VISCOSITY AND THERMAL CONDUCTIVITY TRANSPORT TERMS 

VISCOUS AND THERMAL DIFFUSION TERMS ARE INCLUDED IN THE COMPUTATION OF THE FOREBOOY FLOW AND/OR SHOCKLESS INTERNAL FLOW 
VISCOSITY IS REPRESENTED BY SUTHERLAND'S LAW WITH 

REFERENCE VISCOSITY* 3.500000E-07(LBF-SEC/FT**2I REFERENCE TEMPERATURES 4,920000E402<DEG R) 

BASE TEMPERATURES 1.9e6000E*02(DEG R) 

THERMAL CONDUCTIVITY IS REPRESENTED BY A QUADRATIC CURVE FIT OF TEMPERATURE BASED ON THE FOLLOWING DATA 

DATA POINT NO. 1 TEMPERATURES if , 0 00 QOOE+02 1 DEG R) THERMAL CONDUCTIVITY* 2.550963E-03(FT-SLUG/SEC**3-DEG R> 

DATA POINT NO. 2 TEMPERATURES 1.400000e+03(OEG R) THERMAL CONDUCTIVITY* 7.566417e-0S< FT-SLUG/SEC**3-OEG R> 

DATA POINT NO. 3 TEMPERATURE* 2.400000E+031DEG R> THERMAL CONDUCTIVITY* 1.145772E-02fFT-SLUG/SEC**3-OEG R) 

. Selected output for Sample Case No. 3. 
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Figure 11. Continued. 
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Figure 11. Continued. 


INTERNAL flow field without shock wave ststem 
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5. SAMPLE CASE NO. 4 


Sample Case No. 4 is concerned with the computation of the internal flow 
field in the Mach 3.5 mixed-compression inlet documented in Reference (4). In 
this computation, the internal shock wave system is discretely fitted. This 
sample case considers the calculation of the internal flow field at the design 
point conditions: free-stream Mach number of 3.5, zero centerbody translation, 
and zero angle of attack. 

The data deck listing for this sample case is presented in Figure 12. 

The first card of the data deck is the title card. English units are to be 
used in the computation, hence KUNIT is left equal to its default value of 1 
in namelist LIST!. To perform the internal' flow field integration with dis- 
crete fitting of the internal shock wave system, both KCALL(2) and KCALL(3) 
are retained at their default values of 1 and 0, respectively, whereas 
KCALL(l) = 0 is specified in namelist LIST!. The internal flow field integra- 
tion termination point is x = 6.9 ft, hence XEND(2) = 6.9. Since at zero 
incidence the flow field is axi symmetric, KSYM is left at its default value of 
3. The parameters KVISCY, RCAVG , and KSGL0B in namelist LIST1 are not used in 
this flow field integration option. No restart file operations are to be per- 
formed, hence IPRSTP and KSTART are left at their default values of 0. In this 
computation, however, the shock wave surface normal unit vector components and 
the incident normal Mach number at each shock wave solution point are to be 
output, hence KPRINT = 2 is specified in namelist LIST1 . 

In namelist LIST2, the free-stream Mach number is specified as MFS = 3.5. 
The free-stream pressure and density, specified by PFS and R0FS, respectively, 
retain their default values. The flow is axi symmetric, so both PITCH and YAW 
retain their default values of 0.0. The cowl lip axial station at the design 
point is located at x = 2.86 ft, hence XI = 2.86 is specified in namelist LIST2. 
This inlet has a conical forebody. Consequently, the initial -value plane flow 
property field for the zero angle of attack case is internally generated in the 
computer program. Thus, KIVS and KC0N retain their default values of 1. The 
input parameter ITAPE in namelist LIST2 is not employed in this computation. 

In namelist LIST3, IST0P is left at its default value of 1 since KSYM = 3. 
The specified number of radial stations on both the initial-value plane and 
on each of the internal flow field solution planes is 21, hence JMAXI = 21 and 
JINLET = 21. The input parameter JLIMIT in namelist LIST3 is not used in the 
computation of an internal flow field. 

In namelist LIST4, the- thermodynamic properties GAMMA and R retain their 
default values. No other input parameters in namelist LIST4 are employed in 
this computation 

The contours of the centerbody and the cowl are specified in namelist 
LISTS. For this inlet, the centerbody contour and the cowl contour are de- 
cribed by equation (3) applied to a number of intervals on each contour. 
Consequently, KBASE retains its default value of 0. The number of axial 
stations used for the geometry description of the centerbody and the cowl are 
11 and 14, respectively; hence NCENT = 11 and NC0WL = 14. Thus, the number of 
intervals on the centerbody and the cowl are 10 and 13, respectively. Since 
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SAMPLt Case NO. 4 

SLISTI KCALL<1»=0, XEND(2) s 6.9, KPRINT=2 S 
SLIST2 MFS=3.5» XI=2.86 S 
*LIST3 U«AXI=21. JINLET=21 S 
SLISI4 S 

*U1ST5 NCENTsllt NC0WL=14, KOCENT < 1) =11*1 , KOC.OWL ( 1 ) =14*1 » 

XCENT 1 1 ) = 0.0, 2.798794, 4.0, 4,2, 4,4, 4,55, 4,7, 4 , 9 , 5.5, 6,28, .6,9. 

RCEN7 t 1 ) =11+0.0 , 

AGENT 1 1 ) a 0.0, .493511* ,70532, .7367, .759. .766, .7565, ,7391, 

,6525, .4, 0.0, 

BCENT(1»= .17633, ,17633. .17633, ,144, .052, 0.0, -.0646, -.1295, -.153, 

0 . 0 , 0 ,. 0 » 

CCENTU)= 0.0, 0.0, ,02020035, -.1774997, -.1600005, -.1693327. 

".1615001, -.03499995, .1651873, 0,0, 0.0, 
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6.9, 

KCOWUtl >=14*0.0, 

ACOUL ( 1 ) si ,0 , 1,004168, 1.0051, .9661, .9364, ,9154, .8766, .864, ,8572. ,85, 
.85, .8839, .9227, 0.0, 

BCOWt-<l) s ,01745001, .01745, -.Oil, -.124, -.1942, -,2i3. -.163, -.093, -.0405, 
0.0, 0.0, ,1Q7, .0729, 0.0, 

CC0WK1)= 0,0, -.04926635, -,06500001, -.1664999, -.2859976, .050000,08, 
.3500016, .3049997, ,1075002, 0.0, .1928, .01025012, -.01512487, 0.0, 

PCOWMlJs 0,0, 4.U0418E-03, -.03240740, -.0300005, 1,279984, 
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8.74970E-O3, O.Ot 
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ORIGINAL PAGE IS 
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Figure 12. Data deck for Sample Case No. 4. 
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equation (3) is used to specify the body radius for all intervals,, KDCENT ( I ) = 

1 (1=1 to T1 ) and KDC0WL{I) =1 (I = 1 to 14). Note that except for 

XCENT(ll) and XC0WL{14), the specifications for the other input parameters for 
the last station on either the centerbody or the cowl are immaterial. That 
iss only the XCENT and XC0WL arrays must be specified for NCENT and NC0WL 
elements, respectively. The other arrays require that only (NCENT-1) or 
(NC0WL-1 ) elements be entered. The 11 elements of the XCENT array are thus 

entered. Although not used in the computation, the 11 elements of the RCENT 

array are specified, each element being 0.0. The arrays ACENT, BCENT, CCENT, 
and DCENT are then entered. For each of these arrays, the last (11th) element 
is arbitrarily specified' as 0.0. In a like manner, the arrays XC0WL, RC0WL, 
AC0WL, BC0WL, CC0WL, and DC0WL are entered. -Since this sample case is for 
the design point condition, the centerbody .translation is 0.0 ft; hence, 

D XT RAN in namelist LIST5 retains its default value of 0.0. 

All convergence tolerances and iteration limits specified in namelist 
LIST6 retain their default values. 

No debug output is to be printed, consequently all input parameters in 
namelist LIST7 retain their default values. 

Figure 13 presents selected portions of the output for Sample Case No. 4. 
Sample Case No. 4 required 566 seconds of central processor time on the CDC-6500 
computer (the execution was terminated by a 10,000 line count limit). 
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THE ANALYSIS OF steady three-dimensional flow in SUPERSONIC mixed-compression aircraft inlets 


A9SJHACT 

the flow field in a supersonic mixed-compression aircraft inlet is computed USING the METHOD OF CHARACTERISTICS 
for STEADY THREE-DIMENSIONAL FLOW, THE oOw SHOCK WAVE AND REFLECTED INTERNAL SHOCK WAVE SYSTEMS ARE, COMPUTED 
USING A DISCRETE SHOCK-FITTING PROCEDURE, THE PROGRAM HAS THE CAPABILITY TO INCLUOE THE INFLUENCE OF MOLECULAR 

TRANSPORT ON THE SOLUTION &Y TREATING THESE EFFECTS AS CORRECTION TERMS IN THE CHARACTERISTICS SCHEME. 
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N.A.S.A. GRANT no. N5R-15-Q0S-191 Fi>R THE N.A.S.A. LEWIS RESEARCH CENTER. CLEVELAND. OHIO. THE PRINCIPAL 
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SAMPLE CASE NO. 4 


SPECIF itO COMPUTATION OPTIONS 

1.) INTERNAL, flow FIELD WITH SHOCK WAVE SYSTEF1 


FLOW SYMMETRY 

AXISYMMETRIC flow 


THERMODYNAMIC model 

A YHERMALLY AND CALORICALLY PERFECT GAS IS SPECIFIED WITH 

SPECIFIC HEAT RATIOsI, 40000 GAS CONSTANTS 1 . 71bl&lE+03< FT«LeF/SLUG-U£G R) 


ViSCObITY AND THERMAL COfJUUCJIVlTY TRANSPORT TERMS- 

viscous and thermal diffusion terms are not included in the computation - inviscid and adiabatic flow is assumed 


ORIENIATION AND FREE STREAM DATA 

ORIENTATION - PITCH= 0 . 00000 ( DEGREES ) YAfc= Ol 00 000 IDE'GRtE S > 
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X-VELOCITY= 3;386452E + 03 ( FTXSEC I Y*VELOCITY= 0,* (FT/SECl- Z-9fei-OfcITY= 0 + ' (FT/SEc'l 


Fi gu.re- 1 3.: - S&l ected: output- for' Safnpl e • Case- No 4 



INITIAL VALUE SURFACE 


AN INTERNALLY GENERATED INITIAL VALUE SURFACE IS SPECIFIED AS BEING LOCATED AT X= 2,860000E+00(FT) 

inoex parameters 

ISTOP= I IMAX=40 UMAXI=21 

JINLETC2X 

INTEGRATION termination points 

INTERNAL FLOW FIELD INTEGRATION- TERMINATES AT X= 6.900000E+00(FT) 

centekbodt geometry 


I 

KDCEMT 

XCENT 

rcent 

ACENT 

bcent 

CCENT 

dcent 



(FT) 

(FT) 

(FT > 


(FT**-1> 

<FT**-2) 

1 

1 

0. 

0. 

0. 

..•1-.7£330QE-Ol 

0. 

0. 

2 

1 

2.798794E+00 

0. 

4.935110E-01 

1.763300E-01 

0. 

0. 

5 

1 

4.00000«E+Q0 

0. 

7.053200E-01 

1.763300E-Q1 

2, U2 00350-02 

-3.3675I2E-0I 

4 

1 

4.20000DE+00 

0. 

7.367000E-01 

1.440000E-01 

-1.774997E-01 

-1.750011E-01 

5 

1 

4.400000E+00 

0. 

7.&90000E-01 

5 .2000002-02 

-1.600005E-01 

-5.925896E-02 

6 

.1 

4,S50001)E + 01) 

«. 

7,63Q0U0E-01 

0. 

-1.G93327E-01 

-2,044475E-01 

7 

1 

4 .70000UE+00 

u. 

7.5B50U0E-01 

-6.4&OOOOE-Q2 

-1.G15001E-01 

-2.499656E-03 

S 

1 

4.900D0BE+0U 

o. 

7.391000E-01 

-1.295000E-01 

-3.49999SE-02 

1.712957E-02 

a 

1 

5. 5000002+00 

0. 

t»,525QU0E-0l 

-1.530000E-Q1 

1*6518730-01 

-4, 9236022-01 

10 

1 

6.280000E+00 . 

0. 

4.000000E-01 

0. 

0. 

0. 

11 

1 

&.90000UE+OU 

0. 

0. 

0. 

0. 

0. 


COWL GEOMETRY 


TRANSLATION FROM DESIGN POSITIONS 0. (FT I 


I 

KOCOWL 

XCOWL 

(FT) 

RCOWL 

(FT) 

ACOWL 

(FT » 

* 

3C0WL 

CCOWL 
(FT**-l ) 

DCOWL 

(Ft**'2> 

1 

1 

2.660000E+00 

0. 

i.ouoouue+oo 

1.745001E-02 

0. 

0. 

2 

1 

3.100000E+00 

0. 

1.004186E+00 

1.7450C0E-02 

-4.926635E-02 

4.1104I8E-03 

3 

1 

3.4000002+00 

G. 

1.005100E+00 

-1.100000E-02 

-6.500001E-02 

-3.2407402-02 

4 

1 

4 • 0 OOoOOE. 0 0 

0. 

9.6B1000E-01 

-1.240000E-01 

-1.664999E-01 

-3.0000SDE-02 

5 

1 

4. 2000002+00 

0. 

9.3G4000E-01 

-1.942000E-01 

-2.659976E-01 

1.279984E+00 

6 

1 

4.300000C+00 

0, 

9.154000E-01 

-2.130000E-01 

5.000008E-02 

2.499996E-01 

7 

1 

4.5000002+00 

0. 

8 , 7660 OOE-Ol 

-1 , 630000E -01 

3.5000I8E-01 

-1.251464E-05 

e 

1 

4,600000e+00 

0, 

6.640000E-01 

-9.3QO0OOL-O2 

3,049997E-01 

-5.499974E-01 

9 

1 

4 . 7do OOOE+DO 

0. 

8.572000E-01 

-4.650000E-02 

1.075002E-01 

-7.812535E-02 

10 

1 

5 , 1 0 OOOUE + OO 

0. 

8.500000E-01 

0. 

0. 

0. 

n 

1 

5.600000E+00 

G. 

8.5000 00^-01 

0. 

1 .928000^-01 

-1.144000E-01 

12 

1 

6.10000UE+00 

0. 

6.839000E-01 

1.070000E-01 

1.025012E-02 

-8.612521E-02 

13 

1 

G.500000E+00 

0. 

9.227000E-01 

7.2900D0E-02 

-1.5124S7E-02 

8.749780E-03 

14 

1 

G.90000UE+00 

0. 

0. 

0, 

0. 

0. 


CO 


Figure 13. Continued 
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OF POOR QUALITY 



TT4 


udnveROENCE TOLERANCES. ITER AT ION LIMITS. AND OTHER PARAMETERS 

CONVERGENCE TOLERANCES AND OTHER PARAMETERS 

CRIT* 1)= I.OOqOOOE-OI 
CRXT ( 21= 1.0GOOO&E-07 
CRITI 3)s I.OOOOOOE-OH 
CRITI 4Js 1.000000E-05 
CRITI 5)= l.QOOOODE-04 
CRITt 61= 1.000000E-04 
CRITI 7) = 5.DOOOOOE-OJ 
CRITI 3>s l.OOOOOOE+OO 
CRITI 91 = 1.000000E-04 
CRIT(10)= B,OOOOOOE«01 
CRITllllS 1.000000E-06 
CRIT(12!= 2.000000E-01 
CRITI 13 >s 1 .OOOOD'OE-05 
CRIT ( 1N)= 4.000000'E-OI 
CRIT 1 15> s 1.COOO0OE-04 
CRITUfi’JS l',OOOOOOE-06 

ITERATION limit's 

I TEND! 11=10, 

ITEN0(2)=1O 
ITEND(3)=1'0 
ITENDI4)s20 
ITENOI 51=10 
ITENO(4)=iO 

INPUT SAFETY -FACTORS 9.7500C0E-01 


JEfogrose ;il'3.. ff G,onit;i'n',a:ed;. 



initial data plane x= 2,s60oo(ft) 


I 

J 

V 

2 

M 

0 

P 

HO . 



(FT) 

(FT) 


(FT/SEC) 

( LBF/FT2 > 

(SLUG/FT3) 

1 

1 

.5003 

0.0000 

3.127 

3269.0 

414.22 

5.-3079E-04 

1 

2 

.5293 

0.0000 

• 3.126 

3269.2 

913.65 

5, 3045E-04 

1 

3 

.5544 

0,0000 

3,130 

3269. S 

912.63 

5.2952E-04 

1 

4 

. 579** 

o.oooo 

3.132 

3270.7 

911.29 

5.2810E-04 

1 

b 

. 60F4 

0.0000 

3.136 

3271.6 

909.32 

5.263OE-04 

1 

b 

.6295 

0.0000 

3.139 

3273,2 

907.01 

S.2417E-04 

1 

/ 

.6545 

0.0000 

3.199 

3279.7 

909.39 

5.217&E-04 

1 

e 

.6796 

0,0000 

3.199 

327 o,9 

901.52 

5.1911E-04 

1 

9 

.70«6 

0.0000 

3.159 

3276,2 

396.41 

5.1624E-04 

1 

10 

.7296 

0.0000 

3.159 

3260.2 

395.09 

5.1317E-04 

1 

n 

.7597 

0.0000 

3.166 

3262.3 

391.57 

S.0990E-04 

I 

12 

,7797 

0.0000 

3.172 

3264.5 

367.65 

5.0643E-U4 

I 

13 

,6097 

0.0000 

3.179 

3266,9 

363,93 

5,0277E*0‘» 

1 

IF 

.e290 

o.oooo 

3,166 

32o9 , 9 

379.60 

9, 9890E-04 

1 

15 

. 6596 

0.0000 

3. 199 

3292,0 

375.4J 

4.9479E-04 

1 

16 

.6799 

0.0000 

3.203 

3299,9 

370.79 

4.9042E-04 

1 

17 

,9099 

0.0000 

3.212 

3296.0 

365.69 

9.8573E-09. 

1 

IB 

.9299 

0.0000 

3.222 

3301,3 

360,46 

4.6063E-04 

1 

19 

.9550 

Q.0000 

3,233 

3305.1 

359,58 

9*75010*04 

1 

20 

.9300 

0.0000 

3.296 

3-«9.3 

347.93 

4. 66o2E-04 

1 

U 

1.0050 

o.oooo 

3,262 

3314.5 

340.04 

4.6101E-04 


BASS FLOR BATE for ENTIRE PLANE: 3,S7b50E4-O0tSLirG/SEC) 


Figure 13. Continued 


FQREBOCY FLOW FIELD 


T 

U 

V 

W 

FT 

TT 

(DEG R) 

(FT/SECI 

(Ft/SEC ) 

(FT/SEC) 

(LBF/FT2) 

(DEG R) 

959.7 

3219.3 

567.9 

0,0 

18399,7 

1344,3 

954,6 

322 4,. 0 

541,9 

0,0 

16399,7 

1344.3 

954 , 3 

3228,5 

518.1 

0.0 

16399.7 

1344,3 

453.8 

3232.8 

496.2 

0,0 

18399,7 

1344.3 

453.2 

3237.0 

475.9 

0.0 

16399.7 

1344.3 

952,9 

3241.1 

456.9 

0.0 

16399.7 

1344.3 

451.6 

3245.1 

439,0 

0.0 

16399.7 

1344.3 

950,-7 

3249,1 

422.1 

0.0 

16399.7 

1344.3 

949,7 

3253.0 

4 06.0 

0.0 

18399.7 

1344.3 

446.6 

3256.6 

390.6 

0.0 

18399.7 

1344,3 

447.5 

3260.7 

375,7 

0.0 

18399.7 

1344.3 

44b, i 

3264,6 

361.2 

0,0 

16399, 7 

1344,3 

445.0 

3266,5 

347,1 

0,0 

18399,7 

1344,3 

443. 6 

3272.4 

333,2 

0, 0 

16399.7 

1344,3 

442.1 

3276, 5 

319.5 

0.0 

10399.7 

1344.3 

44Q.6 

3280,7 . 

305.6 

0.0 

18399,7 

1344,3 

458,9 

3265.1 

291.6 

0,0 

16399.7 

1344.3 

457,0 

3289,7 

277,2 

0,0 

18399.7 

1344,3 

435.0 

3294.7 

262.1 . 

0,0 

18399.7 

1344,3 

432,6 

3300.2 

245.7 

0. 0 

16399,7 

1344,3 

429,8 

3306.7 

227.1 

0.0 

18399,7 

1344,3 





Ulto plane at COWL ENTRANCE X= 2. 86000 ( FT I INTERNAL FLOW FIELD WITH SHOCK WAVE SYSTEM 


I 

J 

v 

2 

M 

0 

P 

RO 

T 

U 

V 

W 

PT 

TT 



(FT) 

(FT) 


(FT/SEC ) 

<LBF/'FT2 > 

(SLUG/FTS) 

(DEG 8) 

-(FT/SEC) 

(FT/SEC) 

1FT/SEC > 

(LBF/FT2) 

(DEG ftl 

1 

1 

.5043 

(1,0000 

3.127 

3269.0 

414.22 

5.3079E-04 

454,7 

3219,3 

567,8 

0,-0 

18399.3 

1344.3 

I 

2 

.5318 

8,0000 

3.126 

32t>9 .'3 

413,74 

5,30366-04 

454,6 

3224.4 

539.8 

o.-o 

18400,4 

1344.3 

1 

6 

.5594 

8.0000 

3.130 

3270.0 

412.53 

5,29246- 04 

454.2 

3229-. 3 

514.0 

0.0 

18400.4 

1344.3 

I 

4 

.5889 

0.0000 

3.133 

3271.0 

410.71 

5.27S8E-04 

453,6 

3234.1 

490.4 

D.'O 

18400,4 

1344.3 

I 

5 

.8145 

0.0000 

3,137 

3272.4 

406.41 

5.2546E-04 

452.9 

3238.6 

468,6 

0,0 

18400,4 

1344,3 

1 

b 

• 6420 

0.0006 

3.142 

3273.9 

405,71 

5.2290E-O4 

452.0 

3243.1 

448.3 

0.0 

18400.3 

1344,3 

1 

7 

.6695 

o.oooo 

3.1.47 

3^75.7 

402.70 

5.2020E-04 

451.1 

3247.5 

428,7 

0.0 

18399.7 

1344.3 

1 

8 

.6971 

0.0000 

3.152 

3277.6 

395.37 

5.1712E-04 

450.0 

3251.8 

410.6 

0.0 

16399, 0 

1344.3 

1 

9 

.7246 

0.0000 

3.156 

3279.8 

395.77 

5.1360E-04 

448.8 

3256.1 

393,7 

0.0 

18399,8 

1344,3 

1 

10 

.7522 

0.0000 

3.165 

3262. 1 

391.93 

5.1O23£-04 

447.6 

32b0 ■ 3 

377.3 

0.0 

1839.9,9 

1344.3 

1 

11 

.7797 

o.oooo 

3.172 

3264.5 

367,85 

5.0643E-04 

446.3 

32b4,6 

361,4 

0.0 

16399.9 

1344.3 

1 

12 

.6072 

0.0000 

3.160 

3267.1 

363.52 

5.0259E-04 

444.8 

3260.9 

345,9 

0.0 

18399,6 

1344.3 

1 

15 

.6345 

0.0000 

3.166 

3289.9 

378.94 

4 .98V6E-C4 

443.3 

3273.2 

330.7 

0,0 

16399.8 

1344.3 

I 

14 

.8623 

0,0000 

3,197 

3292,9 

374,05 

4.9349E-04 

441.7 

3277. 7 

315.5 

0.0 

18399.8 

1344.3 

1 

15 

,6698 

o.oooo 

3,208 

329b. 1 

366.83 

4.6856E-04 

439.9 

3262.4 

300,2 

0.0 

18399.7 

1344.3 

1 

16 

a 9174 

0,0000 

3.217 

3299.7 

363,17 

4.8319E-04 

4.38.0 

3287.4 

284.6 

0.0 

18399.6 

1444.3 

1 

17 

,.9449 

0.0000 

3.228 

3303.5 

357.03 

4.773SE-04 

435,6 

3232,6 

266.3 

0,0 

18399,8 

1344,3 

I 

16 

.9725 

0.0000 

3.242 

3306.0 

350.02 

4.7063E-04 

433.4 

3296.5 

250.9 

0.0 

18359.7 

1344.3 

1 

U 

1,0000 

0,0000 

3.259 

3313,5 

341.55 

4.6247E-04 

430.3 

3305.4 

230.8 

0.0 

18399.0 

1344,3 

1 

0 

1,0006 

0.0000 

3.092 

3256.1 

43b. 21 

S.5053E-04 

461.7 

3255.7 

56.6 

-.0 

16371,6 

1344.3 

1 

1? 

1. 0000 

0,0000 

3.092 

3256.1 

436.21 

5.50S3E-04 

461.7 

3255.7 

56,6 

-.0 

18371.6 

1444.3 


SHOCK WAVE POINT PARAMETERS 

I INCIOENT X-COMP. OF Y-COMP. OF i-CQMP. OF 

NORMAL UNIT UNIT UNIT 

MACH NO. NORMAL NORMAL NORMAL 

I X . 112451E+0 0 -2.750988E-01 -9.61H160E-01 -9.746162E-16 

MASS FLOW RATE FC)- ENTIRE PL'ANE= 3 . 82787E+00 ( SLUG/SEC ) 

X-STEP REGULATION PARAMETERS 

LIMITING POINT - 1= 1, J= 1 SAFETY FACTORS 9.T50000E-01 DELTA-X= 7.2571C6E-02I FT > 


Figure 13. .Continued 



Z L L 


SOLUTION PLANE NO 

i. 1 



X= 2.93257 < FT ) 


INTERNAL FLOW 

FIELD WITH SHOCK 

WAVE SYSTEM 



I J 

y 

Z 

M 

C 

P 

RO 

T 

U 

V 

W 

PT 

7T 

ITS 

itl 


tFT 1 

( FT > 


tFT/SEC) 

(LBF/FT2I 

(SLUGRFT3) 

< DEG R1 

(FT/SEC1 

(FT/SEC > 

(FT/SEC 1 

1 LBF/FT2 1 

(DEG R> 



1 1 

.5171 

0.0000 

3,129 

3269,4 

413.40 

5.3D04E-04 

454,5 

3219.7 

567.7 

0.0 

16398.8 

1344,3 

1 

a 

1 2 

,5440 

•.0000 

3.129 

3269.5 

413.36 

5.3000E-04 

454.5 

3224,5 

540.6 

-.0 

16400.8 

1344.3 

1 

2 

1 3 

.5710 

•.0000 

3.131 

3270.1 

412.37 

S.2910E-04 

454,1 

3229,1 

516.2 

-.0 

18400.6 

1344.3 

1 

z 

1 4 

.5980 

..0000 

3.153 

3271.0 

410.75 

. 5.2761L-04 

453.6 

3233,4 

493.4 

-.0 

16400.8 

1344.3 

1 

2 

1 b 

.6250 

.,0000 

3.137 

3272,2 

406.66 

S.2569E-04- 

453,0 

3236.0 

472.3 

-.0 

16400.7 

1344,3 

1 

2 

1 6 

.6521 

-.0000 

3.140 

3273.6 

406.35 

5.2356E-04 

452.2 

3242,2 

452.3 

-.0 

16400,6 

1344.3, 

1 

2 

1 7 

.6792 

.,0000 

3.145 

3275,3 

403.44 

5. 2089E-04 

451.3 

3246,4 

433.7 

-.0 

16400.3 

1344,3 

1 

2 

l a 

.7063 

-.0000 

3.151 

3277.1 

400.25 

5.1794C-04 

450.3 

3250.6 

416.0 

-.0 

16400.0 

1344.3 

1 

Z 

1 9 

.7335 

-.0000 

3.156 

3279.1 

596.90 

5.1464C-04 

449.2 

3254.7 

399.3 

-.0 

16400,0 

• 1344.3 

1 

Z 

1 10 

.7606 

>.0000 

3,163 

3281,2 

393,31 

b.llblE-04 

448.0 

3256.8 

383,3 

-.0 

16400,0 

1344,3 

1 

z 

1 u 

.7078 

-.0000 

3.169 

3233.5 

369,50 

b.0797£-04 

446.8 

3262.8 

367.9 

-.0 

16400.0 

1344,3 

1 

2 

1 12 

, 6150 

.,0000 

3.176 

3265.9 

385.47 

5.0420E-04 

445.5 

3266.9 

352.6 

-.0 

18399,9 

1344,3 

1 

2 

1 13 

.6422 

-.0000 

3.184 

3288,5 

361.20 

5.0020E-04 

444.1 

3271.1 

336,1 

-.0 

18399.9 

1344.3 

1 

2 

1 14 

.6694 

-.0000 

3.192 

3291,3 

376.67 

4.9595E-04 

442.5 

3275.4 

323.5 

-.0 

18399.8 

1344.3 

1 

2 

1 15 

• 8966 

-.0000 

3.201 

3294,3 

371.66 

4.9142E-Q4 

440,9 

3279.7 

306.9 

-.0 

16399,6 

1344,3 

1 

2 

1 lb 

,9238 

-.0000 

3.210 

3297.4 

.366,74 

4.8658E-04 

439.2 

3264,3 

294,0 

- , 0 

16399, 5 

1344,3 

1 

5 

1 17 

.9510 

-.0000 

3,220 

3300,9 

361,16 

4.61AOE-04 

437.3 

3269.1 

279.0 

-.0 

18399,6 

1344.3 

1 

3 

1 IB 

,9761 

-.oooo 

3.232 

3304.7 

355.17 

4.75&6E-04 

435.2 

3294,2 

263,3 

- . 0 

18399, 3 

1344.3 

1 

3 

1 U 

.9792 

-.0000 

3.232 

3504.9 

354.6s 

' 4 . 75ZBE-0 4 

435.1 

3294.4 

262.6 

-.0 

16399,2 

1344,3 

1 

3 

1 D 

,9792 

-.0000 

3.037 

3236.1 

472.67 

5.8287E-04 

472,5 

3235,6 

57.0 

-.0 

18355 • 1 

1344,3 

1 

1 

1 19 

1.0013 

-.0000 

3.063 

3245.6 

46b. 4i 

5.6773E-04 

467,4 

3245.1 

56.6 

-.0 

16371,6 

1344,3 

1 

3 

SHOCK 

WAVE POINT 

parameters 



» 









i 

INCIDENT 

X-COMP, OF 

Y-COMP. 

OF 2 

-COMP. OF 










normal 

UNIT 

UNIT 


UNIT 










MACH NO. 

NORMAL 

NORMAL 


NORMAL 









1 1 

. 133364E+QC 

I -2 i 750968E-01 

-9. 6141606-01 7. 

3263S6E.11 









MASS FLOW RATE FOR ENTIRE 

PLANE= 

3. 83100E +00 (SLUG/SEC) 









COURANT NUMBER = 

..97500 













X-STtP 

REGULATION PARAMETERS 












LIMITING POINT - 

1= I. J= 

21 

SAFETY 

factor= 

9.750000E-01 


DELTA-X= 

6.491595E- 

02(FT) 






Figure 13. Continued. 
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CO 
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Figure 13. Continued 
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Figure 13. Continued. 


ORIGINAL PAGE IS 
• Op POOR QUALITY - 



ro 


INTERPLANAR RESULTS - REFLECTION .WlThj COWL 


INCIDtNT WAVE UPSTREAM AND DOWNSTREAM SHOCK POINTS 


I 

J 

X 

(FT) 

Y 

(FT) 

Z 

(FT) 

M 

a 

(FT/SEC) 

P 

( LBF/FT2) 

RO 

(SLUG/FT3) 

. T 

(DEG fl) 

U 

(FT/SEC) 

V 

(FT/SEC) 

W 

(FT/SECI 

PT, 

( L'BF/FTgi 

TT 

(DEG R) 

1 

u 

4.4181 

.8913 

, OOOO 

2,133 

2767,3 

1849.21 

1.5448E-03 

697,5 

2738.0 

-522.2 

.'0 

16374.5 

1344,2 

I 

D 

4,4101 

,8913 

, oooo 

1,941 

2634,5 

2503.12 

1.9486E-03 

766.5 

2624.0 

-234.5 

.0 

18309.5 

1344.2 


INCIDENT WAVC SHOCK. POINTS 


I INCIDENT 
NORMAL 
MACH NO, 


X-COMP. OF 
UNIT 
NORMAL 


Y-COMP. OF 
UNIT 
NORMAL 


Z-COMP. OF 
UNIT 
NORMAL 


1 


1.15S651E+00 


-3.661369E-01 9.297716E-01 


2.4927S6E-06 


INCIULNT WAVE UPSTREAM AND DOWNSTREAM DODY STREAMLINE POINTS 


1 

J 

X 

Y 

Z 

M 

0 

P 

RO . 

T . 

. u. 

. V 

W 

PT 

TT 



(FT) 

(FT) 

(FT) 


(FT/SEC) 

(L0F/FT2) 

(SLUG/FT3) 

(DEG R) 

(FT/SEC) 

(FT/SEC) 

(FT/SEC I 

(LOF/FT2) 

COCG R> 

1 

u 

4,4181 

,6913 

.oooo 

2.153 

2787.3 

1849,21 

1.544SE-03 

697,5 

2738,0 

-522.2 

.0 

18374,5 

1344.2 

X 

0 

4t4181 

,8913 

.0000 

1.941 

2634.5 

2563,12 

1.9406E-03 

766.5 

2624.0 

-234.5 

.0 

18309.5 

1344.2 


KLFLLL'TEO -WAVE UPSTREAM AND DOWNSTREAM BODY STREAMLINE POINTS 


l 


X 

(FT) 

Y 

(FT) 

Z 

(FT) 

M 

0 

(FT/SEC) 

P 

(L8F/FT2) 

RO 

(SLUG/FT3 J 

T 

(DEG R‘) 

U 

(FT/SEC) 

V 

(F'T/SEC) 

U 

(FT/SEC) 

PT 

(LBF/FT2) 

TT . 
(DEG R) 

1 

u 

4.4101 

,8913 

, oooo 

1,941 

2634.5 

2563.12 

1.9486E-03 

766.5 

2624,0' 

-234.5 

.0 

18309,5 

1344.2 

1 

0 

4 .4101 

.8913 

, oooo 

1.741 

2460.4 

3476.24 

2.4214E-03 

837.0 

2424.7 

-462,4 

' .0 

18256. S 

1344.2 


RtFLECTEO WAVE SHOCK POINTS 


I 


I 


INCIDENT 
NORMAL 
MACH NO, 


X-COMR. OF 
UNIT 
NORMAL 


Y-COMP. OF 
UNIT 
NORMAL 


2-COMP. OF 
UNIT 
NORMAL 


1.142B16E+00 


-6.SB2802E-01’ 


-7.527730E-01 -2.025662E-06 


^rqure;- T3v. 


Gontinue’d'. 



SZL 


SOLUTION PLANE. NO 


55 


X= A • 42525 (FT ) 


internal flow field with shock wave system 


1 

J 

Y 

z 

M 

0 

P 

KO 



(FT) 

(FT > 


(FT/SEC) 

(LBF/FT2 1 

' (SLUG/FT3) 

1 

1 

.7602 

-.0000 

1.746 

2472,3 

3430,69 

2.3946E-03 

1 

2 

.7656 

-.cooo 

1.735 

246b . 4 

34b0 . 77 

2. 4090 E-03 

1 

i 

.7712 

-.0000 

1.736 

2463.9 

3472.91 

2.4145E-03 

1 

4 

.7770 

-.0000 

1.738 

24oS,2 

3464.16 

2.4097E-03 

1 

5 

.7831 

-,oooo 

1.741 

2467.9 

3448.7? 

2.4018E-D3 

1 

6 

■ 7694 

-.oooo 

1.737 

2465.1 

3460.43 

2.41U4E-03 

i 

7 

.7 = 58 

-.0000 

1.736 

24 oU * 5 

5472.2b 

2 • 41 31E-03 

1 

a 

.6026 

-.0000 

1.754 

2479.9 

3364.14 

2 * 3692E-03 

X 

9 

.8095 

.0000 

1.771 

2495.2 

3297.06 

2. 32 61 E-03 

X 

1 U 

.6167 

.oouo 

1. 768 

2492,7 

33l5.4 6 

2.3350E-O3 

1 

u 

.6?40 

.0000 

1 . 765 

2409,9 

3334.75 

2.3448E-03 

1 

12 

.6316 

. GCO 0 

1.788 

2509.9 

5222. is 

2.2386E-03 

1 

14 

. 6394 

. 0000 

1.616 

2533.7 

3090.78 

2 . 22 20E-03 

1 

14 

.8474 

.0000 

1.822 

2530. 5 

30fc4, 63 

2.2069E-03 

1 

lb 

* 6555 

.0000 

1. 620 

2537.4 

3076,77 

2.2134E-0J 

1 

16. 

.6638 

, 0000 

1.838 

2652.8 

2993.99 

2.1732E-03 

1 

17 

.6724 

.0000 

1.676 

2565.1 

2819.8s 

2.0634£-03 

l 

16 

.6814 

.0000 

1.915 

2bl4,J 

2665.14 

2.0030E-03 

X 

U 

,ee5i 

.0000 

1.925 

2621. b 

262 b . 83 

1.9633E-03 

1 

D 

.6651 

.ooco 

1. 723 

2452,6 

3570,75 

2.46b2£-03 

I 

iy 

.6900 

. 0000 

1.729 

2438.5 

3537,01 

2.4509E-03 

SHOCK 

WAVE POINT 

parameters 




1 


INCIDENT 

X-COMP 

. OF 

T-COMP. 

OF 2 

;-COMP. OF 



NORMAL 

UNIT 

UNIT 


UNIT 



PACH NO. 

NORMAL 

Normal 


NORMAL 

1 

i 

•143292E+00 

-6.S&2802E-01 

-7.S27730E-01 -2. 

02149&E-06 


flASS f-LOW Kate FOR ENTIRE PLANEs 3.83503L+00|SLUG/StC) 

CCUHAIMT Nll.%ERs .97600 
X-STtH REGULATION PARAMETERS 

LiMiiiNo point - i= i, j= i Safety factors 9.75000oe-oi 


T 

(OEG R) 

U 

(FT/SEC) 

V 

1 FT/SEC) 

W 

(FT/SEC) 

FT 

(LBF/FT2) 

TT 

(OEG 'R> 

ITG 

itl 

034.8 

2470.0 

108,2 

-.0 

iel46,8 

1343.6 

1 

3 

037.1 

2465.4 

69.9 

-.0 

18123.3 

1343. s 

1 

3 

838,1 

2463.6 

34.6 

-.0 

18109,2 

1343,5 

1 

3 

837.7 

2465.2 

1.3 

-.0 

18103.0 

1343.6 

1 

4 

836,7 

2467.8 

-27,5 

-.0 

18102,1 

1343,7 

1 

3 

838,0 

2464.8 

-41.7 

-.0 

18103,8 

1343.9 

1 

3 

818.4 

2463.9 

-52.0 

-.0 

18106.6 

1344.0 

1 

3 

832, 3 

2476.5 

-81.9 

-.0 

18118.7 

1344,2 

1 

3 

826.1 

2492.9 

-107.0 

-.0 

18130.6 

1344.4 

1 

2 

827.4 

2490.6 

- 102.6 

-. u 

18141,3 

1344.6 

1 

2 

628.7 

2488.1 

-95.3 

-.0 

16153.2 

1344,6 

1 

3 

620.5 

2507.0 

-119.7 

-.0 

16169.3 

1344.9 

1 

3 

810.5 

2529.3 

-148.3 

-.0 

18188,0 

1344.9 

1 

2 

808, -4 

2534.7 

-147.2 

-.0 

18204 • 5 

1345.0 

1 

2 

609. 2 

2535,9 

-133,7 

-.0 

16220.9 

1345,2 

' 1 

3 

802,0 

2548.7 

-145.1 

-.0 

16238.6 

1345.3 

1 

3 

788.7 

2578.7 

-181.4 

.0 

16262.3 

1345.0 

1 

3 

775,6 

2605,8 

-211.4 

.0 

18267,2 

1344.6 

1 

3 

772,3 

2612,7 

-215,9 

. 0 

18295,9 

1344.5 

1 

3 

843.7 

2412.0 

-445.4 

.0 

16242.5 

1344.5 

1 

1 

841.1 

2415.8 

-455.9 

.0 

10256,5 

1344.2 

1 

3 


DELTA-XS 7,251737E-03(FT> 


figure 13. Continued. 


ORIGINAL PAGE IS 
OF POOR QUALITY 



126 


SOLUTION 

PLANE 

NO. 245 



X= 5. 32102 ( FT ) 

I 

J 

Y 

2 

p 

9 

P 

KO 



(FT) 

(FT) 


(FT/SEC) 

(L6F/F.T2I 

(SLUS/FT3) 

1 

1 

.6797 

-.0000 

2.091 

2743.6 

1969.33 

1.6179E-03 

1 

2 

,6076 

-.0000 

2, 039 

2742.2 

1995.33 

1 • 62140-03 

l 

6 

,6960 

-.0000 

2,037 

2740.3 

2002.72 

1.6257Et03 

1 

4 

.7046 

-.0000 

2.054 

2736 , b 

2010.40 

1 .6302E-0? 

1 

5 

• 7 134 

-.0000 

2.092 

273fc . 7 

2016.24 

1 < 634 6E* 03 

1 

6 

,722b 

-.0000 

2.079 

2735,0 

2026 . 0 c 

1 .6395E-03 

1 

7 

,7319 

coco 

2.077 

2733.4 

2035.66 

1.64 39E- 0 3 

1 

8 

,7414 

-.0000 

2.075 

2732.0 

2040.20 

1.6479E-03 

1 

U 

.7421 

-.0000 

2.075 

273l.g 

2040.53 

1.6461E-03 

1 

D 

,7U?1 

-.0000 

2.032 

2700.9 

21F2.79 

1 .7294Q-03 

1 

9 

.7506 

-.0000 

2.022 

2b94. 0 

2216.11 

1.7465E-03 

1 

10 

.7603 

-.0000 

2.016 

2690.7 

2233.23 

1 1 756264-0 3 

1 

11 

,7699 

-.0000 

2.014 

2667,6 

2249., 26 

1 . 7675 C 703 

1 

12 

,7796 

-.oooo 

2.012 

2666.3 

2234. 



Figure 13, .Continued,, 


INTERNAL FLOW FIELD WITH SHOCK WAVE SYSTEM 


T 

(OEG R) 

. U 

(FT/ SEC) 

V 

(FT/SEC) 

V 

(FT/'SE.C) 

PT 

(LBF/FT.2) 

TT 

(DEO R) 

ITG 

ITL 

716.5 

2713.5 

-406.6 

.0 

17944. &' 

1343.2 

1 

3 

717.1 

2715.6 

-331,0 

.0 

17939. B 

1343.0 

' 1 

3 

717.8 

2717.3 

-354. S 

.6 

17936.7 

1342.9 

1 

3 

7ia;6 

2716. a 

-327.7 

.0 

17935.3 

1342.8 

1 

3 

719.3 

2720.2 

-300.2 

*. 0 

17935.7 

1342.3 

1 

3 

720.1 

2721,4 

-272.1 

-.0 

17937.6 

1342,8 

1 

3 

720.8 

2722,5 

-243.6 

-.0 

17941.1 

1342.7 

1 

3 

721.4 

2723.6 

'-214*4 

-.'0 

17946,3 

1342.7 

1 

3 

721.4 

2723.7 

-212.3 

-.0 

17946.7 

1342.7 

1 

2 

735.5 

2667.9 

-265.0 

-.0 

17946,2 

1342.7 

1 

4 

730.6 

2683. 1 

-292,7 

*, 0 

17953.6 

1342.7 

1 

3 

740.1 

2681,6 

-217.9 

-.0 

17962,6 

1342.8 

1 

3 

741,5 

2660.5 

-196,0 

-.0 

17973.8 

1342.8 

1 

3 



6. SAMPLE CASE NO. ’5 
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Sample Case No. 5 is concerned with the computation of the internal 
flow field in the Mach 3.5 inlet discussed in Sample Case No. 4, for the 
off-design conditions of a free-stream Mach number of 2.5, a centerbody 
foward translation of 0.855 ft, and an angle of incidence of 3.0 degrees. 

This sample case also illustrates the use of the program restart option. 

The data deck listings for this sample case are presented in Figures 
14 and 15. Figure 14 presents the data deck used for the initial execution, 
and Figure 15 presents the data deck used for the restart execution. The 
first card in each data deck is the title card. For the initial execution 
data deck, namelist LIST1 is identical to namelist LIST! in Sample Case No. 

4, except for the specified values of KSYM, IPRSTP, and KSTART. Since the 
vehicle is at a nonzero angle of attack, only one plane of flow symmetry 
exists; hence KSYM = 1. The desired number of solution planes to be com- 
puted in the initial execution is 25; hence IPRSTP = 25. Thus, the execu- 
tion will be terminated after plane No. 25 has been printed. Since this is 
the initial execution and a program restart is to ensue, KSTART = 1 is 
specified in namelist LIST1 . This causes the binary restart file to be 
written on TAPE 4 (it is assumed that this file will be stored and then re- 
trieved when the restart execution is performed). 

The free-stream Mach number for this sample case is 2.5; hence MFS = 2.5 
is specified in namelist LIST2. The angle of attack is 3.0 degrees; hence 
PITCH = 3,0. With the prescribed centerbody translation of 0.855 ft, the 
axial location of the cowl lip, and thereby of the initial -value plane, is 
x - 3.715 ft. Note that the origin of the coordinate system remains at the 
forebody tip, hence a forward centerbody translation corresponds to a rear- 
ward cowl translation. Consequently, XI = 3.715 is specified in namelist 
LIST2. For this execution, the results of Jones (3) are used for the initial 
data since the forebody is conical. Hence, the initial -value plane flow pro- 
perty field is externally specified. Consequently, KIVS = 0 is specified in 
namelist LIST2. The bow shock wave is conical, consequently KC0N is retained 
at its default value of 1. The initial -value plane will be entered by a for- 
matted read of the input file; hence ITAPE is left at its default value of 5. 

All input parameters in namelist LIST3 are retained at their default 
values, except for IST0P. As in Sample Case No. 2, 21 circumferential stations 
are specified in the computed sector (half-plane); hence IST0P = 21. 

All parameters in namelist LIST4, which specifies the thermodynamic model, 
retain their default values. 

All parameters in namelist LIST5 are identical to those in Sample Case No. 
4, except for DXTRAN, which specifies the centerbody (cowl) translation. In 
this case, the centerbody has been translated forward by 0.855 ft, or, alter- 
natively, the cowl has been translated backward by 0.855 ft. Hence, DXTRAN = 
0.855 is specified in namelist LIST5. 

All convergence tolerances and iteration limits are retained at their de- 
fault values. Hence no parameters are input in namelist LIST6. 
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6AHPLE CASE NO. 5 

SLIST1 KCAU_U>=0« XENDI2)=6.9, KSYM=l, KPRINT=2 » IRRSTP=25. KST,ART=1 * 
*LIS1 2 MFS=2.5, PITCH=3.0, XI=3.715, KIVS=0 $ 

SL1ST3 IS10PS21 $ 

*LIST4 $ 

•*LIST5 NCENT=11, NCOWL=14, KDCENT ( 1 > = 11*1 , kDCOWL < 1 ) =14*1 » 

XCENT ( 1 ) = 0.0, 2.798794, 4.0, 4.2, 4.4, 4.55, 4.7, 4.9, 5.5, 6.28, 6. : 9, 
HCENI < 1 ) =11*0 .0 » 

ACENT C i ) = 0.0, ,4935U» ,70532, .7387, .759, .763, ,7595, .7 39l> 

• 6325,- .4, 0.0, 

BCENT ( 1 ) = .17633, .17633, .17633, .144* .052, 0.0, -.0646, -,i295‘, -.153, 

0 . 0 , 0 . 0 , 

CCENiT < 1 ) = 0.0, 0.0, ,02020035, -.1774997, -.1600005, -.1693327, 

-.1615001, -.03499995, ,1651873, 0.0, 0,0, 

DCENT'I 1 ) = 0.0, 0.0, -.3367512, -.1750011, -.05925696-, -.2044475, 
-2.499656E-03, ,01712957, -.4923902, 0.0, 0.0, 

XC0WL<1)= 2.66, 3.1, 3.4, 4.0, 4.2, 4.3, 4.5, 4.6, 4.7, 5.1, 5/6, 6.1, 6.'5, 
6.9, 


KCOWUl )=14*0 .0, 

AC OWL ( 1 ) =1 , 0 , 1.004168, 1.0051, .9681, .9364, .9154, .9769, ,.864» .8572, .95, 
.85, .9839, .9227, 0.0, 

BCQWL 1 1 ) s .01745001, .01745. -.Oil, -.124, -.1942, -.2l3, -.163, -i.093, -.0485, 
0,0, Q..0, ,107, ,0729, 0.0, 

CCOWL ( 1 > = 0.0, -.04926635, -.06500001, -.1664999, -.2859976, . 05000008* 
.350UQ18, .3049997, .1075002, 0.0, .1928, .01025012, -.01512497, 0.0* 

OCOWL ( 1 ) = 0.0, 4 . 110418E-03, -.03240740, -.0300005, 1.279964, 

,2499998, -1 . 251464E-05, -.5499974, -.07812535, 0.0, -.1144, -i08&l2521, 
B.74978E-03, 0.0, 

OXTKANsO.855 $ 


SL1S16 $ 

*L1ST7 S 

6, 550547 3333193E -01-1, 5137 063432436E-08 
-9. 5925328260 476E-06 3.1020272481190E+02 
7 . 8361540748 196F- 01-8 « 774 3 8953 30 69 7E-0 7 
-4,1667a757O8460E-O4 3, 093311 1360508E+02 
9. 1 21760 8163243c -01-1. 021 39240 56373E-0b 
-3.044626576O914E-O4 3 . 071974&716077E+02 
1.0407367557829E+00-1.165345&583677E-06 
-3, 5721134419 579E -04 3 . 0427B59240714E+02 
1, 1692974299334E+ 00-1 ,30929931 0B980E- 06 
-.3.33502O9596933E-O4 3 , 0081671 5274 0 0E + 02 
1.29 785b 1040838E+00-1 .4532527634284E-06 
-3.1251496135275E-04 2.96952900 08655E+02 
1,42641 87702343E+OO-1.5972O62139588E- 06 
-2. 93069650261 12E-04 2 .9269256882197E+02 


2.31087410490B8E+03 
4.3206146580673E-04 
2.3181055968066E+03 
4.3134643529953E-04 
2. 32470 90952 54 9E +03 
4.2922U2486969E-04 
2.3311025783340E+03 
4 » 263 0 3 667 3691 bE -04 
2, 33738 1136741 IE +03 
4.22835635260U7E-04 
2.34 3624510 3795E+03 
4. 1094770 3967&2E -04 
2.3499049147152E+03 
4.146474B842468E-04 


4 ,074694 574 39 99E/ 02 
3, 72036066160 03E+02 
3.43032b520‘6562E+02 
3.1866591519535E+02 
2 ,9754472 7952 15E+ 02 
2.767176512438iE+02 
2.6134263226085E+02 


-1*41 189996673 07E +00 

0. 

-1. 50650562090 55E+ 00 

0. 

-1.60U112750804E+00 

0. 


0. 

3.298152O370OO6E+O2 

0. 

3, 1977663586671E+02 

0. 

3.0670442486558E+Q2 


2.3259929l6l670E+03-8,4506243644906E+oi 

4.5i4o391336456E-04 

2 • 33645048B4367E+03-5 , 78970 869B7370E+ 01 
4.41S47041&5349E-04 

2. 3490464328 799E+ 03-2; 5922170 d48908E+oi 
4.28577191 ;2761E-04 


Figure 14 


Data deck for Sample Case No. 5, initial execution. 
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SAMPLE CASE NO. 5 - CONTINUED 

SUISU KCALL<1>=0i XENO<2)=6,9, KSYM=1, KPRINT==2« IPRSTP=50» KSTART=2 * 
•LIST 2 MFS=2.5, PITCH=3.Q, *1=3.715, KlySsQ $ 

•11ST3 IST0P=21 * 

seism $ 

•LiS»5 NCENT-11, NCOWLs=l<f, KOCENT ( 1 ) =11*1 , kDCOWL< 1 ) =14*1 « 

XCENT(l)s O.O, 2,798794, 4,0, 4.2, 4,4, *.55, 4.7, 4.9, 5,5. 8. 28, 6.9, 
RCENl (1)311*0,0, 

ACLNT ( 1 ) = 0.0, ,493511, ,70532, .7387, ,759, ,763, .7585, ,7391, 

•6525, ,4, 0.0, 

BCENKDe ,17633, .17633, .17633, .144, .052, 0.0, -.06*6, -.1295, -.153, 
O.o, 0,0, 

CCENT ( 1)3 0.0, 0,0, .02020035, -.177*997. -,1600005* -.1693327, 

“.1615001, -.03499995, ,1651873, 0,0, 0.0, 

DC£Nltl)= 0,0, 0,0, -,3367512, -.1750011, -,05925696, -.20*4*75, 
-2.499656E-03, ,01712957, -.4923802, 0,0, 0,0, 

XCOWL(l)= 2.86, 3.1, 3.4, 4.0, 4.2, *.3, 4.5, 4.6, *.7, 5.1, 5.6, 6.1, 6, 
6,9, 

KC0WU1 ) =14*0 , 0 , 

ACOWL(1)=1.0, 1.004188, 1.0051, ,9681, .936*, ,9154* ,8768, ,864, .8572, 
.85, .8839, .9227, 0.0, 

BCOWL <11— .01745001, .01745, -.011, -.124, -,i942, -,2l3* -.163, -,093, - 
U.O, 0,Q, .107, .0729, 0.0, 

CCOWL < 1 ) = 0,0, -.04926635, -.06500001, -.1664999, -, 2859976, .05000008, 
.3500018, ,3099997, .1075002, 0.0, .1928* .01025012, -,01512487, 0.0, 
DCOWUlla 0.0, 4, 110418E-03 « -.03240740, -.0300005, 1,279984, 

.2499998, -1 .251464E-05, -.5499974, -.07812535, 0,0, -.1144, -.08812321, 
6.74978E-03, 0.0* 

UXTKAN=o.655 i 
$L1ST6 $ 

SLI6T7 $ 


5, 

.89 , 
.0485, 


Figure 15. Data deck for Sample Case No. 5, restart execution. 



No debug output is to be printed. Hence, all parameters in namelist 
LIST7 are left equal to their default values. 

The data deck for the initial -value plane flow property field follows 
namelist LIST7. The complete listing of the data deck for the initial -value 
plane is not shown to conserve space. Since the bow shock wave is conical 
(the forebody is conical), the BETA angles are internally generated and are' 
not specified by the user (K.C0N = 1). 

Figure 16 presents selected portions of the output for the initial exe- 
cution of Sample Case No. 5. This execution required 755 seconds of central 
processor time on the CDC-6500 computer. At the end of this execution-, the 
restart file written on TAPE 4 was stored. 

Figure 15 presents the data deck- for the restart execution. This data 
deck is identical to the initial execution data deck, except that the initial - 
value plane flow property field is not entered, and the specified values for • 
IPRSTP and KSTART in namelist LIST! are different. The initial -value plane 
for this execution is solution plane No. 25 from the initial execution., which 
is read- from the' binary restart file written on TAPE 4 (TAPE 4 is linked- to 
the dummy file RESTRT in the PR0GRAM card). The new value for the' plane 
number at which the program execution is to be terminated is 50; hence 
IPRSTP' = 50 is specified in namelist LIST!. Since this is a program restart,- 
KSTART = 2 is specified in namelist LIST! . All other namelists are identical 
to those used in the initial executi on o f Sample Case No. 5. 


Figure 17 presents selected portions of the output for the restart exe- 
cution of Sample Case No. 5. The solution terminated jus-t after solution 
plane No. 46. This execution required 590 seconds of central processor time 
on the CDC-6500 computer. 
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THE ANALYSIS OF STEADY THREE-DIMENSIONAL FLOW IN SUPERSONIC MIXED -COMPRESSION AIRCRAFT INLETS 


ABSTRACT 


THE FLOW FIELD IN A SUPERSONIC M I XEd-C OmPRESSI ON AIRCRAFT INLET IS COMPUTED USING THE METHOD OF CHARACTERISTICS 
FOR STEADY THREE-DIMENSIONAL FLOW. the BOW shock wave ano reflected internal shock wave systems are computed 

* DISCRETE SHOCK ^FITTING PROCEDURE* THE PROGRAM HAS THE CAPABILITY TO INCLUDE THl INFLUENCE OF MOLECULAR 

TRANSPORT ON THE SOLUTION bY TREATING THESE EFFECTS AS CORRECTION TERMS IN THL CHARACTERISTICS SCHEME, 


THIS PROGRAM WAS DEVELOPED AT THE PURDUE UNIVEHSITT THERMAL SCIENCES AND PROPULSION CENTCR 
N.A.S.A. GRANT NO. NfaR - 1 5 - 00 S>- 19 1 FoR THE N.A.S.A. LEWIS RESEARCH LENTER. CLEVELAND t OHIO. 

investigator was j.d, Huffman and the n.a.n.a. technical director was a. bishop. 


BY J. VADYAK UNDER 
THE PRINCIPAL 


JOB TITLE 

SAMPLE case NO, 5 


SPECIFIED COMPUTATION OPTIONS' 

1,1 INTERNAL flow field with shock wave system 


FLOW SYMMETRY 

ONE PLANE OF SYMMETRY - COMPUTED SECTOR IS THE HALF-PLANE BOUNDED BY THE Y-AXIS AND CONTAINING THE +Z-AXIS 


thermodynamic mooel 

A THERMALLY and calorically perfect gas is specified with 

SPECIFIC hCAT RATIO=l, 40000 GAS CONSTANT= 1.71G161E + OMFT-LBF/SLUG-OLG R) 

VISCOSITY AND THERMAL CONDUCTIVITY TRANSPORT TERMS 

VISCOUS AND THERMAL DIFFUSION TERMS ARE NOT INCLUDED IN THE COMPUTATION - INVISCID AND ADIABATIC FLOW- IS ASSUMED 


okieniation ano free stream data 

YAW= 0.00000<DEGREES) 

PRESSURES 2,422000E+02{LBF/FT**21 DENSJTY= 3, 622C0GE-D4 I SLUG/FT**3 I 

TEMPERATUftEc 3.896437E+02 IDEG R! SONIC SPEEU= 9. 675S77E+02 (FT/SEC ) 

X-VELOCITYs 2.415579E-t-03(FTKSEC! Y-VELOCITY= 1 .265951E+02 ( FT/SEC 1 2-VELOCITY= D. (FT/SECI 


ORIENTATION - PITCHs 3.00000 (DEGREES! 
FREE STREAM DATA - MACH NO.= 2.50000 


Figure 16. Selected output for Sample Case No. 5, initial execution 
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INITIAL VALUE SURFACE 

A USER SUPPLIED INITIAL VALUE SURFACE IS READ FROM TAPE NO. 5 AND IS SPECIFIED AS BEING LOCATED AT X= 3. 715000E+00 ( FT) 


index parameters 

ISTOP=21 IHAX=40 JMAXI=11 

UIKLETzII 


integration termination points 

internal FLOW FIELD INTEGRATION TERMINATES AT X = 6.900000E+00 (FT) 


CENTEKfaODT GEOMETRY 
I KDCENT 


1 1 

a i 

i i 

4 1 

b I 

6 1 

7 1 

6 1 

9 1 

10 1 

11 1 


XCENT RCEtJT 

(FT)' (FT) 

0 . 0 . 

2.798794E+0Q 0. 

4.00DQUOE+OU 0, 
4.200000E+00 0. 

4 , 400000^400 0 . 

4.5b0000e+00 0, 

4.70000UH4.0U 0, 

4.900000E+00 0. 

5.bU00O0£+OO 0, 

fe.aaooooE+oo o, 

6,90 OOOOE 4 OO 0, 


ACENT 

(FT) 

0. 

4.93511OE-01 
7. 0S321UIE-01 
7, 3870 llOE-Ul 
7.SSU0UOC-01 
7,b30000£.01 
7,b85U0G£.Ol 
7.391000E-01 
6,b2S0yOE-01 
4,tlOOOOUE»Oi 

a. 


SCENT 


J.76330OE-01 

I.763300E-01 

1.763300E-01 

1,4AOOOO£401 

b.200000£'02 

0. 

-6.460000E-02 

-1.29b00O£-0l 

-J..5300UDE-U1 

0. 

0, 


CCENT 
( FT**-1 > 

0, 

0. 

2.020035E-02 
-1.774997E-01 
-1 , 60000SE-01 
-1.693327E-01 
•I.falSOOlL-Ol 
-3.499995E-C2 
1.651873E-01 
0. 

0. 


CCENT 

<FT**»2> 

a. 

0 . 

-3.387512E-01 

-1.750011E.01 

-5.925696E-02 

-2.044475E-01 

-2.499656E-03 

1.7129b7E-02 

-4.923B02E-D1 

0. 

0. 


cowl geometry 


TRANSLATION FROM DESIGN POSITIONS 8 . 550000E-01 t FT 1 


I 

kocowl 

XCOWL 

(FT) 

RCOWL 

(FT) 

ACOWL 

fFT) 

SCOWL 

CCOWL 

1FT**-1) 

DcOWL 
( FT**-2) 

1 

1 

2.06OOOOE+OU 

0. 

1, OOOOUUB+OO 

1.745001E-02 

0, 

0, 

2 

1 

3.100000E4DO 

0. 

l.ooRiaBE+oo 

1 .745000E-02 

-4,92663bE-02 

4.11041BE-03 

3 

1 

3.400000E.+ 00 

0. 

1. 005100E+00 

-1 , 100000L-02 

-6.500001E-02 

-3.240740E-02 

4 

1 

4.000000E400 

0. 

9.68100CE-01 

-1.240000E-01 

-1.664 999E-0i 

-3.000050E-02 

5 

1 

4.2UOOUUE+OU 

0. 

9.364000E-01 

-1.9420UOE-U1 

-2 . 859976E-01 

1.279904E+OO 

6 

1 

4.300000E+00 

u. 

9.1540UOE-01 

-2,l30U0OEW01 

b,OOOOOBE-02 

2.499998E-01 

7 

1 

4 . bOOOUVE+OU 

u. 

8, 7680 WOE- 01 

-1 , 630000E-01 

i.b00018£-01 

-i,25i«f64E-05 

8 

1 

4 ,600'0C0L»00 

0. 

8, 8400 0 0E-01 

L9.300000E-02 

3.04*997E-6i 

-S.-4999'?9E-C1 

9 

1 

4.700000C+00 

0. 

6.S72000E-U1 

-4*'850000E-02 

ij-uTSdoeEioi 

-tiSiESASE-’OE 

10 

1 

S.IOOOOUE+CO 

0 . 

8, bOOOUUE-Ol 

6.' 

Ui . 

0. 

11 

1 

5.6OO0OUE>OU 

0 . 

8, SOoOOOE-tll 

0 . 

i.gasoootAoi 

- 1 . 144000E-01 

12 

1 

6.100000E+00 

0 . 

8, 839000E-01 

1 , 07OOOQE-O1 

l,O2S0i2e-0'2 

- 8 , 612 S 21 EL 02 

13 

1 

6 . bOOQOOE+O0 

0 . 

9, 22700 OEWtfi 

7.290000C-02 

F1J512487E-02 

8.74970OE-O3 

14 

1 

6.90000U£4t>0 

0 . 

0 . 

0, 

0 , 

0; 


Figure 16. Continued. 



CONVENIENCE TOLERANCES i ITERATION LIMITS , ANO OTHER PARAMETERS 


CONVERGENCE TOLERANCES ANO OTHER PARAMETERS 


CRITC 11 = 
CRITI 21= 
CRIT ( 31= 
CRIT ( 41= 
CRIT £ 51 = 
CRITI 61= 
CRITC 71= 
CRITC a)= 
CRITI 91= 
CRIT(101= 
CPlTtll 1 = 
CRIT (12(= 
CRIT ( 13 I = 
CRI T C 14 1 = 
CRIT ( 15 1 = 
CRIT l 16 1= 


l.OODOOOE-Ol 
1.000000E-07 
1.000000E-D4 
l.OOOOOOfc-OS) 
I.OOOOOOE-C4 
1.000000E-04 
5.000000E-01 
1 • 00000 OE+OO 
l.OOOUOOE-OR 
6.000000E-01 
l.OQOOOOE-Ofc 
2.000000E-01 
l.OOCOOOE-05 
4 • OOOOOOE” 01 
1 ( DOOOOOE-04 
1.0UOOOOE-Q6 


ITERATION limits 

itendui=io 
I TENO I 2 1 =10 
I TENO (31=10 
ITEND (41=20 
ITENU (51=10 
ITENO(61=10 

INPUT SAFETY FACTOR= 9.750000E-01 


_ Figure 16. Continued. 
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rer 


INITIAL DATA PLANE 


X: 3*71500 ( FT ) 


FOREBODY FLOW FIELD 


1 

J 

r 

2 


a 

P 

■KO 

T 

U 

V 

W 

PT 

TT 



(Ft) 

IFT) 


(FT/SEC) 

(LBF/FT2* 

(SLUG/FT3) 

(DEG R) 

( FT/SEC > 

(FT/, SEC), 

(FT/SEC) 

.(LBF/FT2) 

t DEO R) 

1 

1 

.6551 

-.0000 

2.341 

234b. S 

310.20 

4.32O6E-04 

418.4 

2310.9 

407,. 5 

-.0 

.4132.6 

87,6.7 

2 

1 

.6970 

.1025 

2,340 

2346.3 

310.42 

4.3227E-04 

418.4 

2310.4 

407.5 

32.0 

41-32', 8 

076,7 

3 

1 

.6230 

.2024 

2.339 

2145,6 

311.10 

4.3296E-04 

418,7 

2309.2 

406.4 

.66,9 

4132,7 

,076:7 

4 

a 

• 5837 

,2974 

2.33b 

2344.5 

312,27 

4 .341 1£ -04 

419.1 

2307,1 

4.0.3,, a 

104,0 

4132,8 

876.7 ' 

5 

i 

,5300 

.3850 

2.333 

£342,8 

313.97 

4.3560E-04 

419,8 

2304.3 

397,6 

144,0 

4132,6 

■876.7 

6 

i 

, 9632 

,4632 

2.328 

2340,6 

3lb,28 

4^38092-04 

420.7 

2300.9 

386.0 

186 . 0 

4132.7 

876.7 

7 

i 

.3850 

.5300 

2.322 

2337,7 

319.20 

4.40982-04 

421.6 

2296.8 

367.0 

233 \ 9 

4132,6 

876.-7 

a 

i 

.2979 

.5837 

2.315 

2334,2 

322.62 

4.44552-04 

423.1 

2292.3 

336.6 

261.2 

4132. .6 

876.7 

9 

i 

.2029 

.6230 

2.307 

233U.1 

327.08 

4,48722-04 

424.7 

226,7.3 

300.3 

32b. 5 

4132.6 

876.7 

lu 
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Figure 16. Continued. 
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Figure 16. Continued. 
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REDISTRIBUTED PLANE AT CO-L ENTRANCE Xs 3 , 71500 <FT> INTERNAL FLOW FIELD WITH SHOCK WAVE STSTEM 
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190.6 

4133,2 

6 76,7 

10 

u 

- ,8910 

,4340 

2.245 

2299.3 

560,10 

4 .8071C-04 

456.6 

2280,7 

-163,0 

149.6 

4152.7 

876 .7 

19 

u 

-.9511 

.3090 

2.241 

2297.4 

362.29 

4.8273E-04 

437.3 

2287,5 

-107,1 

103.2 

4132.6 

876.7 

20 

u 

-.9877 

.1664 

2,239 

2296.2 

363,57 

4.6394E-04 

437,0 

2286.6 

-201.8 

52.5 

4132,3 

876.7 

21 

u 

-1,0000 

• 0 0 0 0 

2.2.78 

2295,9 

364.01 

4.6436E-04 

437,9 

2266,5 

-207,0 

.0 

4132.4 

676,7 

1 

0 

1,0000 

o.oooo 

2.077 

2208.4 

4b5,7 1 

5.7652E-04 

470,7 

2208,1 

38,5 

-.0 

4106,2 

876,7 

2 

0 

,9877 

.1364 

2.077 

22O8.5 

465,69 

5.7&49E-04 

470,7 

2206.0 

41,6 

-1&.4 

4106,4 

67b,7 

3 

u 

.9511 

,3090 

2.077 

2208.8 

4bS,3& 

5.7623L-04 

470.6 

2208,0 

50.5 

-50.7 

4107.1 

876.7 

It 

0 

.6910 

,4540 

2.078 

2209.1 

465.02 

5, 7597E-04 

470.4 

2207,6 

64.1 

-40,9 

4106,2 

676.7 

b 

D 

.8090 

,3876 

2,079 

2209.7 

464,45 

3.7532E-04 

470.2 

2207.8 

80,9 

-45,8 

4109,6 

876,7 


U 

.7071 

,7071 

2,060 

2210.2 

4 o4 .04 

5.7523E-04 

470.1 

2207.5 

90.0 

-44.3 

4111.3 

876.7 

7 

0 

.5678 

.8090 

2.061 

2213.6 

463,49 

5.7401E-04 

469.0 

2207.4 

115.6 

-36.5 

4113.1 

876.7 

a 

u 

.4540 

.8910 

2.061 

2211.1 

463 , 29 

5.74 71E-04 

469,7 

2207.2 

129,6 

-22,8 

4115,1 

876,7 

9 


,30 90 

,9511 

2.062 

2211.5 

4e2.99 

5,74526-04 

4b9.6 

2207.2 

130.9 

-4. 6 

4117.0 

876.7 

10 

u 

,1564 

,9677 

2,062 

2211,6 

463,14 

5.74 7 3L-04 

469,6 

2207, U 

141,6 

16.6 

4118.0 

876,7 

11 

u 

.oooo 

1,0000 

2 , 0o2 

2211 , b 

463.23 

5,74876-04 

469.5 

2207, U 

138.2 

38.5 

4120,4 

876,7 

12 

D 

-.156“ 

,9677 

2 . Doc 

2211 .4 

463. 7« 

3 . 7S34£-g4 

469.6 

2206,9 

127.3 

59,2 

4121.8 

676.7 

13 

0 

-.3090 

,9511 

2.061 

2211.2 

464.13 

5,7S/7£-b4 

469,. 7 

2207.1 

111.0 

7b. b 

4123.1 

67b. 7 

14 

0 

-,45«0 

,6910 

2,061 

2210-0 

464 ,7b 

5.7638E-U4 

469,9 

2207.2 

89.2 

88. 7 

4124,0 

376.7 

lb 

D 

-.3678 

.6090 

2,060 

2210.4 

465.39 

5.769SE-04 

470.0 

2207,4 

65.0 

94,8 

4124,9 

876,7 

lb 

□ 

-.7071 

.7071 

2.075 

2210,0 

465.9b 

5.7747E-C4 

470.2 

2207,6 

39,0 

93.5 

4125,3 

676.7 

17 

D 

- .6090 

.3876 

2.0&0 

22ib.i 

4 ob. 56 

5.7715E-04 

4 70,0 

2206.5 

14,2 

85.2 

4123,4 

S7b , 7 

16 

D 

-.8910 

.4540 

2.079 

2210,0 

465,91 

5.7747E-04 

470.1 ' 

2206,8 

-7,4 

70.3 

4125,6 

076.7 

19 

U 

- .9511 

,3090 

2,079 

2209,7 

466,34 

5.770b=-U4 

470,2 

2208,9 

-29.2 

50.5 

4125,0 

676,7 

20 

D 

-.9877 

,1564 

2.079 

2209,6 

46b. 42 

5.7793E-04 

470.3 

2209.1 

-34.9 

2b, 1 

4125,8 

876,7 

21 

0 

-1 ,0000 

.oooo 

2.076 

2209.4 

466.62 

5.7011E-04 

470,3 

2209.1 

.38.5 

.0 

4126,0 

676.7 

1 

9 

l.cnoo 

0.0000 

2.077 

2206,4 

465,73 

5.7652E-04 

470.7 

2206.1 

38.5 

-.0 

4106,2 

876,7 

2 

9 

,9877 

,1564 

2.077 

2206.5 

465,69 

5.7649E-04 

47.0,7 

2206.0 

41.6 

-16,4 

4106,4 

076,7 

3 

9 

.9511 

..3090 

2.077 

2203,6 

465,36 

5.7623E-C4 

470,6 

2206,0 

50,5 

-30.7 

,4107,1 

876,7 

4 

9 

.6910 

.4540 

2 » 076 

2239.1 

465,02 

5.75970-0,4 

470.4 

2207.6 

,64.,1 

-40,9 

4108,2 

876,7 

b 

9 

.6090 

.3B78 

2.079 

2209,7 

464,45 

5.7552E-04 

4.70.2 

2207.8 

80, 9 

-45,8 

4109,6 

076,7 

6 

9 

.7071 

.7071 

2.080 

2210.2 

464,04 

5.7523E-04 

470.,1 

2207.5 

98. .6 

-44.3 

4,111,3 

076,7 

7 

9 

.5870 

,8090 

2.081 

2210.8 

463,49 

5.7481E-04 

469.0 

2207.4 

115.8 

-36.5 

4113.1 

07,6,7 

& 

9 

..4540 

,8910 

2.081 

2211. .1 

463,29 

5.7471E-04 

469.7 

2207.2 

129.6 

-22*8 

4115,1 

876,7 

9 

9 

.3090 

,9511 

2.962 

2211.5 

462,95 

5.7452C-04 

469.6 

2207.2 

130.9 

-4.6 

4117, ,0 

676.7 

10 

9 

.1564 

.9B77 

2.082 

2211.6 

463,14 

5 ,74 73E-.04 

469.6 

220.7.0 

1.4,1, .6 

,16.6 

4118,8 

076.7 
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HtOISTRigUTED 1 

=lane at 

COWL ENTRANCE 

X= 3 # 71500 (FT) 

I 

j 

Y 

z 

K 

0 

P 

’ KO 



(FT) 

(FT) 


(FT/SEC) 

(LBF/FTa ) 

(SLUG/FT 3 ) 

11 

9 

,0000 

1,0000 

2.082 

2211.6 

463. 23 

5 . 7487 E -04 

12 

9 

-.1564 

.9677 

2,082 

2211.4 

463.70 

5 , 7534^-04 

15 

9 

-.3090 

.9511 

2.081 

2211 .2 

484 . 1 $ 

5 . 7577 E -04 

l 4 

9 

-i 4540 

.8910 

2.081 

2210.8 

464,76 

5 . 7636 E -04 

lb 

9 

-.5878 

.8090 

2.060 

2210.4 

465.39 

5 » 7695 E - 04 

16 

9 

-.7071 

.7071 

2.079 

2210.0 

4 b / 5,96 

O. 7747 E -04 

17 

9 

-.8090 

,5878 

2.060 

2210.2 

465.56 

5 . 7715 E -04 

Ifl 

9 

-.8910 

.4540 

2.079 

2210.0 

465,91 

5 . 7747 E -04 

1 ? 

9 

-.9511 

.3090 

2.079 

2209.7 

466,34 

5 . 7766 E -04 

2 U 

9 

-.9877 

. 15 b 4 

2 . 079 

2209.6 

466,42 

5 . 7793 E -04 

21 

9 

- 1.0000 

.0000 

2.076 

2209.4 

46 o.b 2 

3 . 78116-04 


INTERNAL flow field with SHOCK WAVE system 


T 

II 

V 

W 

PT 

TT 

(DEG R) 

{FT/SEC ) 

(FT/SEC ) 

(FT/SEC ) 

(LBF/FT 2 ) 

( OEG R) 

469,5 

2207,0 

138,2 

38,5 

4120,4 

876,7 

469.6 

2206.9 

127.3 

59,2 

. 4121.8 

876.7 

469.7 

2207.1 

111.0 

76.6 

4123.1 

676.7 

469.9 

2207,2 

89,2 

88,7 

4124.0 

876.7 

470,0 

2207.4 

63.0 

94,8 

4124.9 

876.7 

470.2 

2207,6 . 

39,0 

93,5 

4125,3 

676,7 

470.0 

2206.5 

14,2 

85.2 

4125.4 

876,7 

470.1 

2208.8 

- 7.4 

70.3 

4125.6 

876,7 

470,2 

2208; 9 

. 24,2 

50.3 

4125,8 

876,7 

470,3 

2209.1 

- 34.9 

26,1 

4125.8 

876.7 

470,3 

2209.1 

- 36.5 

.0 

4126,0 

676,7 


shock wave point parameters 


l ’INCIDENT 
normal 

‘ MACH flO. 


X-COhP, OF 
UNIT 

normal 

-3.8a2S22£’-01 

-3.88392BE-01 

-3.8B6944t-01 

-3.697054L-Q1 

13.9UB3B7E-01 

-3.922443t-oi 

-3. 9392355^01 

-3.93B213E-01 
-3.979l64i'-0l 
-4.00lA14E-vl 
-4.Oa4399t-.Cl 
-4.047624L-01 
-4.0704686-01 
-4 . 092591E-0 1 
-4 i 11262 OL- 01 
-4.130680L-01 
-4 • 14 70266-0 1 
-4.149669E-C1 
-4 ,lb877bt-01 
-4.174&05E-01 
-4 • 1 762Bbt-Q i 


T-COMP. OF 
UNIT 
NORMAL 

-9.215616E-01 
-9.ioi4aaE-oi 
-6,7bl9l6E-01 
-8,205&i2E-01 
-7 • 44 o671£-01 
-6.SU4400E-01 
-5.402S76C-01 
-4. 169120LF01 
-2,8349891-01 
-1 , 453656E -01 
-1 . 27 43U9L-14 
1.440471E-01 
2. 6225B1E -01 
4.142533E-01 
5.457764E-01 
6.4395S9L-01 
7.3617U7E-01 
8.102633E-01 
8.644749L-Q1 
8.97bl21E-01 
9.086179L-O1 


2 -COMP. OF’ 
UNIT 
NORMAL 

- 9 . 3440 & 3 E- 1 & 
- 1 . 441534 E -01 
- 2 . 046919 E-U 1 
- 4 .l 8097 OE. 0 I 
-S. 410323 E -01 
- 6 . 50440 UE -01 
- 7 . 436 QC 8 E -01 
* 8 . 1823 S 8 E - 01 
- 8 . 725199 E -01 
- 9 . US 1747 E -01 
- 9 . lb 4464 E -01 
- 9 . 031636 L -01 
- 8 . 667011 E -01 
- 8 . 129766 L -01 
” 7 < 374329 L-Q 1 
- 6. 4 39 S 59 E- U 1 
- 5 * 34 Qb 9 $E -01 
- 4 . 128498 E -01 
- 2 . 808849 E-U 1 
- 1 . 421519 E -01 
- 7 . 93288 /E-lb 


1 1 « 204 35PE + O0 

2 1.203829E+00 

3 1.2020016+00 

4 ‘ 1.199209L+00 
b 1.193219E+00 
b 1.190522E+00 

7 1,1848261+00 

8 1.178773L+00 

9 1 . 172 012L+00 
1U l.lb&246E+00 

11 1 1 1581471+00 

12 1 . 151 3C 9'1 + 00 

13 1.1446074+00 

14 1.138358E+00 
lb 1.152645L+00 

16 1 ,127bl‘3t + U0 

17 1 . 122 310L + 00 

18 1 .U 8752 E + 00 

19 1 . 116329L+00 

20 1 . 114 bisL +00 

21 1 * 114 279L + 00 


MASS FLOW RATE FOR ENTIRE PLANCr 1 .86292E+0 0 ( SLUS/SEC 1 
X-STEP REGULATION PARAMETERS 
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Figure 16. Continued. 
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X= 3.78810(FT> 


internal flow field with shock wave system 


X 

J 

Y 

Z 

M 

a 

P 

KO 

T 

U 

V 

W 

PT 

TT 

ITG 



(FT) 

(FT) 

• 

(FT/SEC) 

(LSF/FT2) 

(SLUG/FT3) 

(DEG R) 

(FT/SCC ) 

(FT/SCC! 

(FT/StC) 

(L6F/FT2) 

(DEG R) 


1 

1 

.btro 

-.0000 

2.341 

239 6, <3 

309,65 

4.3171E-04 

413.2 

2311,2 

4 07» 5 

-.0 

4132.8 

876,7 

1 

2 

1 

.6599 

.1035 

2.391 

239b, 7 

310.04 

4.31V0E-04 

418,3 

2510,8 

407.3 

32.6 

4132,0 

876,7 

1 

A 

1 

.6357 

,209a 

2.390 

2345,1 

310.68 

4.0253E-U4 

418,5 

2309.7 

406.4 

66.7 

4132,8 

876.7 

1 

4 

1 

.5965 

.3007 

2.357 

2395.0 

311.74 

4.3359E-C4 

416,9 

2307,7 

4 03,6 

103,3 

4132,7 

876,7 

1 

S 

1 

. 5426 

.3696 

2.339 

2393.4 

313.33 

4.3517E-09 

419. b 

2305,0 

397,5 

143.2 

4132,6 

876,7 

1 

b 

1 

.4755 

.9692 

2,330 

2341.3 

315.47 

4.3727E-04 

420.4 

2301,7 

386.1 

186,5 

4132.5 

876,7 

1 

r 

1 

.39*7 

, 5379 

2.329 

2336.6 

318.24 

4.90U3E-U4 

421.4 

2297,7 

367,4 

232,4 

4132,4 

.076,7 

1 

b 

1 

.H 6 « 

,5926 

2.317 

2335,3 

321,68 

4.4342E-04 

422.7 

2295,5 

339.7 

279,2 

4132,3 

876,7 

1 

y 

1 

.2120 

,6339 

2.309 

2331.3 

325.77 

4.4749E-04 

424.2 

2268,6 

301.6 

324,6 

4132.2 

876.7 

1 

10 

1 

. 110 b 

, 6567 

2.500 

2326 . 6 

330.52 

4.5209E-O4 

426,0 

2263,9 

252.8 

365.9 

4132.1 

8 / 6 , 7 

1 

ii 

1 

, 0061 

,667V 

2 . 

2321.7 

335,61 

4.5724E-U4 

427, -V 

2278,8 

193.3 

400.1 

4132.0 

076.7 

1 

la 

1 

096b 

,6606 

2.279 

2316,3 

391,56 

4.6202E-O4 

430.0 

2273,6 

124,7 

424,0 

4131.9 

876.7 

1 

l* 

1 

-.2CP7 

,6370 

2.260 

2310, 6 

347,55 

4, 68*lt-04 

432,2 

2268. 7 

48.9 

434.9 

4131. 8 

676. 7 

1 

14 

1 

- » 2 9ch 

, 5976 

2.257 

2o03. 2 

353.66 

4. 7447E-G4 

434.3 

2264.3 

-30.7 

430.9 

4131.7 

876.6 

1 

lb 

1 

3eo7 

,6932 

2,296 

2297, 8 

359.5b 

4, B012L-04 

436,4 

2260 . 1 

-110.9 

419.6 

4131.7 

876, 6 

1 

lb 

1 

- , 4hV3 

.9753 

2.236 

2279.6 

3bb.l2 

4.864IE-04 

438.3 

2256,'3 

-187.6 

373.9 

4131.7 

076.6 

1 

1 / 

1 

-.6363 

,3959 

2.228 

2270.9 

370.01 

4, 7004E-04 

4«Q,Q 

2253,1 

-257,0 

321,2 

4131,6 

87b, 6 

1 

16 

1 

-,593S 

,3056 

2.221 

2266.7 

374.11 

4, 9392E-04 

441,4 

2250.5 

-315,4 

254.6 

4131,6 

87b. 6 

1 

iy 

1 

-.634 7 

.2061 

2.216 

22»4.l 

377. 13 

4.7678E-04 

4 42,8 

2246.7 

-359.4 

176,5 

4131,6 

676.6 

1 

2u 

1 

-.£596 

,1059 

2.212 

2262,4 

379.05 

4, 9B57E-04 

443.0 

2247,5 

-306,9 

90.3 

4131,6 

076.6 

1 

21 

1 

- . b6£0 

,0000 

2.211 

22tl, fl 

375.65 

4, 9714E-04 

443,2 

2247.2 

-376,2 

.0 

4131,6 

076.6 

1 

1 

2 

,710? 

-.0000 

2.391 

239b, b 

310.14 

4.3206E-04 

418.3 

2315.2 

395.3 

-.0 

4134,9 

876,7 

1 

2 

2 

.7021 

.1102 

2.391 

204o.5 

310.36 

4.3231E-04 

416.4 

2312.8 

394.9 

32.0 

4135,0 

876,7 

1 

6 

2 

. £>765 

.2179 

2.539 

fc046.e 

011.09 

4 , 3301E- 04 

*tlh t 6 

2311.7 

393,4 

65.2 

4135.0 

376,7 

1 

4 

2 

,63u4 

,3202 

2,337 

2394.6 

312.32 

9. 3423E-04 

419*1 

2309,8 

390.0 

100.6 

4135.0 

676.7 

1 

b 

2 

.5770 

,919<). 

2.533 

2o4 2 , 9 

314. 07 

4. 3577L-04 

419,8 

2307,2 

303.1 

13B.9 

4134,9 

876,7 

1 

b 

2 

.5055 

,9999 

2.326 

2340.6 

314.93 

4,38305-09 

420,7 

2304.0 

371.1 

lao.o 

4134,9 

876,7 

1 

7 

2 

.9216 

,5720 

2,322 

2307.7 

319.40 

4.4123E-04 

421 .8 

2300,2 

352,2 

223,1 

4134.8 

876.7 

1 

d 

2 

,3?73 

,qS06 

2.315 

2304.2 

323.03 

4,94811-04 

423.2 

2296.0 

324,9 

866.0 

4134.7 

876,7 

1 

5 

2 

.2299 • 

,6739 

2,307 

2330.1 

327.29 

4,46 78E-04 

424,6 

2291,5 

200.0 

408,7 

4134,6 

676,7 

1 

1U 

2 

, 110.9 

,7007 

2.297 

232a, 5 

332,15 

4 . 5373E-04 

426,6 

2266,8 

241.1 

346,7 

4134,5 

876,7 

1 

11 

2 

.nOSS 

.7105 

2,267 

2320.4 

337.51 

4.5895E-64 

426,5 

226? • 0 

164. 6 

377,6 

4134.3 

676,7 

1 

12 

2 

1059 

,7029 

2,276 

2315,0 

343.27 

4,64 52E-U4 

430.6 

2277,2 

120.0 

398,9 

4134,2 

876.7 

1 

Lb 

2 

-.2192 

.6772 

2.265 

2307.5 

34y,24 

4,70271-04 

432,7 

2272.6 

49,2 

408.0 

4134.1 

876.7 

1 

14 

2 

-.3179 

,6351 

2.259 

2304 , 0 

355.23 

9, 7601E-04 

434.8 

2268.3 

-25,3 

403.2 

4i33,9 

676. 7 

1 

IS . 

, 2 

-.9137 

,5772 

2.299 

2298.7 

361.05 

4.81560-0“ 

436.9 

2269.3 

-99.7 

383,5 

4133.6 

876. 7 

1 

16 

2 


.5099 

2.234 

2c93.0 

368.43 

4 , 86876-04 

4 38, T 

2260.8 

-170.6 

348.4 

“133,6 

876.7 

1 

1/ 

2 

-.5725 

,9201 

2.226, 

2289, b 

371.22 

4.9120E-04 

440.9 

2257.8 

-234.9 

299.0 

4133.5 

076.7 

1 

16 

2 

-.6315 

,32“ b 

2.215 

2256.1 

575.14 
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Figure 1 6„ Continued: 
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Figure 16. Continued. 
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-9. 10 2387E-0 1 

-1.4 S5891E-01 

3 

1 , 209522L + 0 0 

-3.88&876E-01 

.8,7b5416E-01 

-2.8J6215E-01 

4 

1.20681SC.OO 

•3 , 896906L.01 

-8.21 3133E-0 1 

-4.166S62E-01 

a 

1.203101E+00 

-3.90ei36E-01 

-7.459158E-01 

•5.393277E-01 

6 

1.198889E+00 

-3.922077E-O1 

-6. 522265C-0 1 

-6.48670BE-01 

7 

1. 193459L+00 

-3. 9388016*01 

-5.4Z3&82E-01 
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8 

1 , 187916E .00 
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1.126627E+00 
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8. 9744000-01 

-1.426114E-01 

21 

1.128503E+00 
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9.0861790*01 

6.82345'1E-11 


MASS FLOW RATE FOR ENTIRE PLANE= 1.B8362E+00 ( SLUG/SEC ) 
COUKANT NUMBERS .97500 


X-STtR REGULATION PARAMETERS 

LIMITING POINT - 1= 5, J*U SAFETY FACTOR: 9.7S0'0O0E-Dl 
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16 
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-.5047 
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17 
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-.5983 
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1.973 

2147,9 

546.58 
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INTERNAL flow field with shock wave system 
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INTERNAL FLOW FIELO WITH SHOCK WAVE SYSTEM 
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THE ANALYSIS OF STEADY THREE-DIMENSIONAL FLOW IN SUPERSONIC MIXED-COMPRESSION AIRCRAFT INLETS 


ABSTRACT 

the FLOW FIELD IN A SUPERSONIC MIXED-COMPRESS: ON AIRCRAFT INLET IS COMPUTED USING THE METHOO OF CHARACTERISTICS 
FOR STEADY THREE-DIMENSIONAL flow, ThE fOw SHOCK WAVE AND REFLECTED INTERNAL SHOCK WAVt SYSTEMS ARE COMPUTED 

USING A DISCRETE SHOCK-FITTING PROCEDURE. THE PROGRAM HAS THE CAPABILITY TO INCLUDE TH(1 INFLUENCE OF MOLECULAR 

TRANSPORT ON THE SOLUTION RY TREATING These effects as CORRECTION TERMS in the CHARACTERISTICS scheme. 

IHIS PROGRAM WAS DEVELOPED AT THE PuRCUE UNIVERSITY THERMAL SCIENCES AND PROPULSION CENTER 8* J. VAOYaK UNDER 
N.A.S.A. GRANT NO. NGR-1S-00S-191 FOR THE N.A.S.A. LEWIS RESEARCH CENTER. CLEVELANO. OHIO, THE PRINCIPAL 

INVESTIGATOR WAS J.O. HOFFMAN AND THE N.A.S.A. TECHNICAL DIRECTOR WAS A. BISHOP. 


JOB TITLE 

SAMPLE CASE NO, 5 - CONTINUED 


SPECIFIED COMPUTATION OPTIONS 

1.1 INTEhNAL FLOW FIELD WITH SHOCK WAVE SYSTEM 


FLOW SYMMETRY 

ONE PLANE OF SYMMETRY - COMPUTED SECTOR IS THE HaLF-PLaNE BOUNDED BY THE Y-AXIS AND CONTAINING THE +Z-AXIS 


THERMODYNAMIC MODEL 

A THERMALLY AND CALORICALLY PERFECT GAS IS SPECIFIED WITH 

SPECIFIC HEAT RATI0=1 .40000 GAS CONSTANT:: 1 . 71G161E + 03 ( FT-LBF/'SLUG-UEb R| 


VISCOSITY AND THERMAL CONDUCTIVITY TRANSPORT TERMS 

VISCOUS AND TKEPRAL DIFFUSION Tt-RrlS ARE NOT INCLUDED IN THE COMPUTATION - INVISCID AND ADIABATIC FLOW IS ASSUMED 


OHIENTATION and FREE STREAM DATA 

ORIENTATION- PITCHs 3, 0 0000 ( DEGREES) 
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Figure 17. Selected output for Sample Case No. 5, restart execution 
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